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Approximate controllability via adiabatic 
techniques for the three-inputs controlled 
Schrodinger equation. 

Francesca Carlotta Chittaro * Paolo Mason ' 


Abstract 

We consider a system described by a controlled bilinear Schrodinger equation with three 
external inputs. We provide a constructive method to approximately steer the system from 
a given energy level to a superposition of energy levels corresponding to a given probability 
distribution. The method is based on adiabatic techniques and works if the spectrum of the 
Hamiltonian admits eigenvalue intersections, with respect to variations of the controls, and if 
the latter are conical. We provide sharp estimates of the relation between the error and the 
controllability time, and we show how to improve these estimates by selecting special control 
paths. 


1 Introduction 


A typical issue in quantum control concerns the controllability of the bilinear Schrodinger equation 



Hr, 


Z^ 

fc =i 


u k {t)H k ip(t), 


( 1 ) 


where ip belongs to the Hilbert sphere of a (finite or infinite dimensional) complex separable 
Hilbert space H and H 0 ,... ,H m are self-adjoint operators on H. Here H\. .... H m represent 
the action of external fields on the system, whose strength is given by the scalar-valued controls 
Mi,..., u m , while Hq describes the uncontrolled dynamics of the system. 

The controllability problem aims at establishing whether, for every pair of states ipo and ip\ 
in the Hilbert sphere, there exist controls u k (-) and a time T such that the solution of (jTJ) with 
initial condition ip(0) = ipo satisfies ip(T) = ip\. 

While the case where H is a finite dimensional Hilbert space has been widely understood [41F4], 
in the infinite dimensional case the answer is far from being given. In particular, negative results 
have been proved when TL is infinite-dimensional (see mm)- Hence one has to look for weaker 
controllability properties as, for instance, approximate controllability (see for instance Em OSH 
[25]I. or controllability between subfamilies of states (in particular the eigenstates of Hq, which 
are the most relevant physical states) or in more regular subspaces of square-integrable functions 
(see [6J[7]). 

While the above mentioned works are essentially obtained by means of non-constructive argu¬ 
ments, the purpose of this paper is to propose a method that permits to explicitly select control 
inputs steering the system from the initial state to an arbitrarily small neighborhood of the given 
target state. Adiabatic theory and conical intersections between eigenvalues constitute the main 
tools of the control strategy we propose in this paper. 
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Roughly speaking, the adiabatic theorem (see mmm) states that the occupation probabil¬ 
ities associated with the energy levels of a time-dependent Hamiltonian H(-) are approximately 
preserved along the evolution given by iip(t) = H(t)ip(t), provided that H(-) varies slowly enough. 
This result works whenever the energy levels (i.e. the eigenvalues of #(•)) are isolated for every 
t. On the other hand, if two eigenvalues intersect, and provided that H(-) is smooth enough, the 
passage through the intersections may determine (approximate) exchanges of the corresponding 
occupation probabilities (see [301 Corollary 2.5] and [TO]). For these reasons, adiabatic methods 
are largely used in quantum control to induce population transfers (see for instance the techniques 
known as Stimulated Raman Adiabatic Passage (STIRAP), Stark-chirped rapid adiabatic pas¬ 
sage (SCRAP)) and to prepare superposition states [20]. The applications of adiabatic methods 
in quantum control, as a tool for obtaining controllability results, have already been exploited 
in previous papers (see for instance IZUEBEZlEg). The general idea is to use slowly varying 
controls, taking advantage of the adiabatic theorem, and “climb” the energy levels through the 
conical intersections. 

Related to the present paper is a method recently developed in m in the case m = 2 and for 
self-adjoint Hamiltonians with real matrix elements. It exploits a generalization of ESI Corollary 
2.5] stating that it is possible to arbitrarily recombine the probability weights associated with two 
subsequent energy levels by following (slowly) a suitable control path passing through a conical 
intersection between them. The control strategy of m applies whenever a part of the spectrum 
of the Hamiltonian operator is uniformly separated from the rest of the spectrum (as a function 
of the control parameter), is discrete and each pair of subsequent eigenvalues intersect in a conical 
intersection. When there exists such a portion of the spectrum, called separated discrete spectrum, 
this control strategy permits to attain (approximately) a state having a prescribed distribution 
of probability (relative to the energy levels of the separated discrete spectrum) starting from an 
eigenstate. In particular this entails a controllability property, that we call spread controllability, 
which, although weaker than the usual approximate controllability property, is more practical. 
Note indeed that the relative phases between pairs of components in the eigenbasis decomposition 
are essentially uncontrollable since they evolve according to the gaps between the corresponding 
energy levels. Furthermore, notice that this method allows us to control the population inside 
some portion of the discrete spectrum, if well separated from the rest, even in the presence of 
continuous spectrum, unlike many other classical methods. 

Concerning the precision of the method, an application of the adiabatic theorem together 
with m Corollary 2.5] shows that the maximal error is of the order of the square root of the 
control speed. On the other hand in m it was shown that the precision of the transfer may be 
remarkably improved if one follows some special paths in the space of controls; namely, such paths 
permit to attain a state with a prescribed probability distribution with an error of the order of 
the control speed. From a practical point of view this means that, to guarantee a given precision, 
one may significantly reduce the duration of the process, whose extent constitutes one of the main 
disadvantages of the implementation of adiabatic techniques. 

The purpose of this paper is to adapt the control strategy introduced in nsa to the general 
case of self-adjoint Hamiltonians, assuming that three controlled Hamiltonians are employed, and 
to select the control paths that allow to improve the precision of the process as explained above. 
Preliminary results in this sense were discussed in mi- Notice that the chosen setting is quite 
natural, since it is well-known, for Hermitian matrices or within spaces of self-adjoint operators 
satisfying particular transversality conditions, that the set of operators admitting multiple eigen¬ 
values is a submanifold of codimension three (see e.g. [3, |3T] 33]). Moreover conical intersections 
do not constitute a pathological phenomenon since, as shown in Appendix [B] all eigenvalue in¬ 
tersections are generically conical in the finite dimensional case and in some physically relevant 
infinite dimensional models. Conical intersections are also structurally stable with respect to vari¬ 
ations of the Hamiltonian operator, as shown in Theorem |4.8[ Concerning the relation between 
conical intersections and controllability properties of the bilinear Schrodinger equation, let us fi¬ 
nally mention the main results of the recent paper [12]: if all subsequent energy levels of the 
Hamiltonian are connected by means of conical intersections then the system is approximately 
controllable and, in the finite dimensional case, it is even exactly controllable. Notice that these 
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results have not been obtained by adiabatic techniques, although, as shown in Section [5] it is not 
difficult to recover approximate controllability of the system by extending the results of this paper 
(in a non-constructive way). 

The structure of the paper is the following. In Section [2] we introduce the notations used 
throughout the paper, the main assumptions and definitions, and we adapt the classical statement 
of the adiabatic theorem to our setting. In Section [3] we discuss some properties of conical inter¬ 
sections and related results that allow to propose the basic control strategy. In Section [4] we define 
special paths and, by means of a series of technical results, we show that they can be included in 
the control algorithm in order to improve its performance. As a byproduct, we get a structural 
stability result concerning conical intersections. In Section [5] we briefly mention some extensions 
of the control strategy and of the controllability results obtained earlier. Appendix [A] gathers the 
technical results concerning the regularity of the spectrum and of the spectral projections that are 
needed throughout the paper, while Appendix [B] discusses the genericity of conical intersections 
in the finite and infinite dimensional cases. 


2 Notations and preliminary results 

We start this section by introducing the notations that will be used in the rest of the paper. 

For a function /(•) of a real parameter s, we use the following notation for its right and left limits 
at so: 

/(«o ) = lim , /(«)■ 

s—>S 0 

Moreover we say that f(s) = o(s k ) if lim^o = 0. 

Whenever y(s), s G [si, S 2 ] is a curve on R 3 and Q(-) is a function of v G R 3 then, with abuse of 
notations, we denote by Q(j(r)) the derivative of the composition Q(j(-)) computed at r, that is 

Q(l(r)) ~ f s Q( 70))| s =r = ^(7(r))S0)- Similarly, Q ( %(r)) := |JrQ(7(s))U=r- 
The scalar product of two elements if 1 , 1^2 in the Hilbert state space is denoted by (if 1 , 1 ( 2 ), while 
the scalar product of two vectors wi,w 2 in any other euclidean space is denoted by wi • w 2 . 
Analogously, the norm in the two cases is denoted respectively by || • || and | • |. 

For a given vector v or matrix A the respective transpose is denoted by v T and A T . The inverse 
of the transpose of an invertible square matrix A is denoted with A~ T . 

Given a vector v = (vi,v 2 ,v 3 ) G C 3 , we denote its complex conjugate (v* , v 2 , v 3 ) by v* and its 
real and imaginary parts respectively by 

fKc(v) = (fHe(ui),9le(u 2 ),fHe(u3)) 3m(v) = (3m(ui), 3m(u 2 ), 3m(w 3 )). 

The symbol id is used to denote the identity operator on a vector space which is specified at each 
occurrence, whenever not clear from the context. 


2.1 General setting 

Let H be a separable complex Hilbert space with norm || • |: let us introduce the following notion 
of relative boundedness between operators: 


Definition 2.1 (A-smallness and A-boundedness) Let A, B two densely defined operators with 
domains V(A) C V(B). We say that B is A-bounded if there exist a,b > 0 such that \\Bif\\ < 
a||A0|| + b\\i/>\\ for every if G T>(A). B is said to be A-small if for every a > 0 there exists f3 > 0 
such that ||-B'0|| < allA^H + /3||t/'|| for every ■i/i G T>(A). (The latter notion is called infinitesimal 
smallness with respect to A in m) 


Given a self-adjoint operator A on "H, for every A-bounded operator B we define its norm with 
respect to A as 


\\B\\a 


\\m\ 

WA) Ill’ll + W 


( 2 ) 
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This provides a norm in the space C(T>(A), H) of continuous linear operators from 'D(A) (endowed 
with the graph norm of A) to H. 

We consider the Hamiltonian 

H(u) = H 0 + uiHi + u 2 H 2 + u 3 H 3 , 

with u = (ui,u 2 , u 3 ) € R 3 , and where Hi, i = 0, ..., 3 satisfy the following assumption: 

(HO) Hq is a self-adjoint operator on a separable complex Hilbert space H, and Hi are ffo-small 
self-adjoint operators on H for i = 1,2,3. 

Under assumption (HO), i27; Theorem X.12] guarantees that H( u) is self-adjoint with domain 
V(H 0 ). Moreover, it is easy to see that for every u, H 0 is ff(u)-bounded, and therefore Hi is 
H (u)-small, for every i = 1,2, 3, with constants a, b (as in Definition [2d]) that depend continuously 

on u. 

Schrodinger Hamiltonians are typical Hamiltonian operators describing quantum phenomena 
and can be represented in the form — A + V on the Hilbert space L 2 (H), where D is a domain of 
R n , A is the Laplacian on fl (with Dirichlet or Neumann boundary conditions) and V : D —> R 
has to be interpreted as a multiplicative operator on L 2 (H). In particular such Hamiltonian 
operators are unbounded operators. In this context Hypothesis (HO) is thus intended to describe 
an Hamiltonian operator of the previous form that can be controlled by means of three external 
inputs so that H 0 = — A + Vo and Hi = Vi for some multiplicative operators V), for 0 < i < 3. 

Finite dimensional representations of quantum systems are also common, for instance in the 
description of spin systems. In this case the Hamiltonian operator H( u) is a Hermitian matrix. 
Consider for instance the case of a spin-1/2 particle immersed in a controlled magnetic field. In 
this case, Hi are the Pauli matrices, and the controls are the components of the magnetic field. 

The dynamics of the quantum systems we consider are described by the time-dependent 
Schrodinger equation 

(3) 

Such an equation has mild solutions under hypothesis (HO), u(-) piecewise C 1 and with an initial 
condition in the domain of H 0 (see e.g. [27:, Theorem X.70] and |5j). 

We are interested in controlling <[3j) inside some portion of the discrete spectrum of H( u). Since 
we use adiabatic techniques, some spectral gap condition is needed: 

(HI) There exist a domain u> in R 3 , a map X defined on w that associates with each u £ oj a 
subset E(u) of the discrete spectrum of H( u), and two continuous functions fi, f 2 : w —t R such 
that 

• /i( u ) < /a(u) and E(u) C [/i(u), / 2 (u)] Vu e u. 

• there exists T > 0 such that 


inf inf 

u£o; AESpec(if (u))\£(u) 


dist(A,[/i(u),/ 2 (u)])) 


> r. 


In this case we say that E(u) is a separated discrete spectrum. 

Notation From now on we label the eigenvalues belonging to a separated discrete spectrum E(u) 
in such a way that E(u) = |A 0 (u),..., Afc(u)}, where A 0 (u) < • • • < Afc(u) are counted according 
to their multiplicity (note that the separation of E from the rest of the spectrum guarantees that 
k is constant). Moreover we denote by </>o( u ), - - -, </>fc( u ) an orthonormal family of eigenstates 
corresponding to Ao(u),..., Afc(u). Notice that in this notation Ao does not need to be the ground 
state of the system. 

Our techniques rely on the existence of conical intersections between the eigenvalues, which 
constitute a well studied phenomenon in molecular physics (see for instance mmmmm )- in 
this paper we will adopt the following definition, consistent with the one already given in for 
the two-inputs case (Figure[l]shows a conical intersection in this latter setting). 
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Figure 1: A conical intersection for a controlled Hamiltonian with two inputs. 


Definition 2.2 Let H(-) satisfy hypothesis (HO). We say that u £ R 3 is a conical intersection 
between two subsequent eigenvalues A j and A J+1 if Xj(u) = Aj_)_i(u) has multiplicity two and there 
exists a constant c > 0 such that for any unit vector v £ R 3 and t > 0 small enough we have that 


A J+ i(u + tv) — Aj(u + tv) > ct. 


(4) 


A discussion on the occurrence of conical intersections and on their genericity in some relevant 
cases is provided in Appendix [B] 

To conclude this section, let us make some remarks on the regularity properties and the asymp¬ 
totic behavior of the eigenfamilies of H( u) in our setting. Notice that in general the regularity 
properties of the Hamiltonian induce similar regularity properties of the eigenfamilies, see Propo¬ 
sition A.4 In particular, thanks to the Lipschitz continuity of the eigenvalues, @ holds true in a 
neighborhood of a conical intersection, that is there exists a suitably small neighborhood U of u 
and C > 0 such that 


Aj + i(u) — Aj(u) > C|u — u|, Vu £ U. (5) 

Moreover, it is well known that the eigenvectors can be chosen analytic along straight lines u(-) 
possibly passing through eigenvalues intersections (see [IE],[251 Theorem XII.13]). 

Consider a C 1 curve u : I — > R 3 and assume that the eigenvalues A/ : I —> R and the 
eigenstates : I —» LL, l = 0,..., k are C 1 (/). By direct computations we obtain that for all t £ I 
the following equations hold: 

A i(t) = (4>i(f), (wi(f)#i +u 2 {t)H 2 + u 3 (t)H 3 )$ l (t)} (6) 

(A m {t) - A i{t)) ($;(t), <f> m (f)} = ($;(t), (tn(t)Ffi + ii 2 {t)H 2 + u 3 {t)H 3 ')$ m (t)). (7) 

If Aj(u) = A i+1 (u), then, thanks to |7|, for every half-line r v (t) = u + tv with v = (ui,U 2 ,U 3 ) 
unit vector and t > 0, we have 


^hm + (<hj(r v (f)), (viHi + v 2 H 2 + v 3 H 3 )$ j+1 (r v (t))) = 0. 


( 8 ) 


2.2 The adiabatic theorem 

In this section we recall a classical formulation of the time-adiabatic theorem ([El [Hi M M) 
adapted to our framework. For a general overview see the monograph [13- 

Let H( u) = H 0 + Yli=i u iHi satisfy (HO)-(Hl). Assume that the map / = [to,t/] 9 t >->■ 
u(t) = (ui(t),u 2 {t),u 3 {t)) belongs to C 2 (/,K 3 ). We introduce a small parameter e > 0 that 
controls the time scale, and the slow Hamiltonian H(u(et)), t £ [t 3 / e,Tf/e\. In this notations, 
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r is a geometric parameter used to describe the curve in the space of controls, while t is the 
“chronological” time of the evolution along the control path. 

We denote by U £ {t,to) the time evolution (from t 0 = To/e to t = r/e) generated by H(u(et)), 
and with U £ (t,to) the time evolution generated by the Hamiltonian H a {et), where H a {r) = 
if(u(r)) — *£P*(u(t))P*(u(t)) — *£P^ l (u(t))P*- l (u(t)) is the adiabatic Hamiltonian , P*(u) denotes 
the spectral projection of H( u) on E(u), and Pfi (u) = id — P* (u). 

Theorem 2.3 Assume that H( u) = Hq + satisfies (HO)-(HI). Let I C R and u : 

/ —> u> be a C 2 curve. Then P* £ C 2 (I,C(TL)) and there exists a constant C > 0 such that for all 
tq , r £ I, and setting to = To/e,t = r/e, 

II U e (tfi 0 ) - U £ (t,t 0 )\\ < Ce (1 + e\t - f 0 |). (9) 

Remark 1 If there are more than two parts of the spectrum which are separated by a gap, then it is 
possible to generalize the adiabatic Hamiltonian as (\2fj ) H a {r ) = P(u(r))— ie'}2 a P a (u(T))P Q (u(T)), 
where each P a (u(r)) is the spectral projection associated with the separated portion of the spectrum 
labeled by a. 


Let us now assume that E = {Ay, Aj+i}; we can take advantage of the adiabatic theorem to 
decouple the dynamics associated with the band E from those associated with the rest of the 
spectrum, in order to focus on the former. 

Let W(t) denote the subspace spanned by the eigenstates associated with Aj(u(r)) and 
A_/-|_i(u(t)). Since W(r) is two-dimensional for any r, it is possible to map it isomorphically 
on C 2 and identify an effective Hamiltonian whose evolution is a representation of U £ (t, to)|w(ei 0 ) 
on C 2 . In particular, if we can find a C 1 eigenstate basis {$i(u(t)), $ 2 (u(r))} of W(r) (asso¬ 
ciated with a reordering {Ai(u(r)), A 2 (u(t))} of {Aj(u(r)), A J+ i(u(r))}), then the isomorphism 
U(t) : W(t) —>• C 2 is continuous. Represented in C 2 , the evolution U £ (•, to)lw(et 0 ) i s governed by 
the Hamiltonian Hf s (et), where H £ s (-) is the effective Hamiltonian, whose form is 

H e M = f A i( u W) 0 ^ _ A $ i( u ( t ))-^i( u ( t ))) ( $ 2 (u(r)),$ 1 (u(r)))\ , 

effl J 'v 0 A 2 (u(t)); ^(^(uM),^^))) ($ 2 (u(r)),$ 2 (u(r)))j’ 1 J 


with associated propagator Uf s (t,to) =U(et)U £ {t,to)U*{£to). 
Theorem 2.3 implies the following. 


Theorem 2.4 Assume that {Aj,Aj+i} is a separated discrete spectrum on some oj £ R 3 and let 
u : I —» u> be a C 2 curve such that there exists a C 1 -varying basis o/W(-) made of eigenstates of 
H( u(-)). Then there exists a constant C such that for all tq,t £ I, and setting to = tq/e, t = r/e, 


|| (U e (t.,t 0 ) -U* {et)U £ efi (t,to)U(£to)) | w(et 0 ) II < Cs(l+e\t - t 0 |). 


3 Conical Intersections and general control strategy 

3.1 Properties of conical intersections 


Conical intersections have a characterization in terms of the non-degeneracy of a particular matrix, 
which contains some geometric properties of the eigenspaces relative to the intersecting eigenvalues, 
as shown below. 


Definition 3.1 We define the conicity matrix associated with two orthonormal elements £ 

V(Ho) as 


M(i/j 1 ,^ 2 ) = 


{i’i 1 H 1 if 2 )* 
{ 1 / 1 , H 2 i> 2 )* 
(tfi,H 3 4> 2 )* 


(V’pPiV^) 
(V’l, H 2 lf>2) 
{if\, H 3 if 2 ) 
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(ip 2 , H-l^) - (■j/'ljPl^l) 
( if 2 ,H 2 xf> 2 ) - {'fi-i_,H 2 xf> 1 ) 
{i>2,H 3 1p 2 ) - (V'l,P'3V’l) 



Lemma 1 The quantity det ip 2 ) is purely imaginary and the function (ipi, V’ 2 ) H > det 1 P 2 ) 

is invariant under unitary transformation of the argument, that is if (ipi,ip 2 ) T = U('0i, ip 2 ) T for 
a pair ip i,ip 2 of orthonormal elements of V(Hq) and U € U(2), then one has det ip 2 ) = 


detM(ipi,ip 2 )- 


Proof. The fact that the determinant is purely imaginary comes from direct computations. To 
prove its invariance, we set 


/gd/h+fo) coga gdft+fe) s i nQ A 

\ — e _ ^ 2 sino! e _J/3l cosa J 


( 11 ) 


for some real scalars /3i, /? 2 , (3 3 , so that ip 1 = e 1 ^ 3 (e 1 ^ 1 cos a ipi + e 1 ^ 2 sin aip 2 ) and^ 2 = —e~ 1 ^ 2 sinaipi + 
e _J/Sl cos aip 2 . 

By direct computations it follows that 


/ g- 2 */3i-*/33 cog 2 a 

M(tpi,ip 2 ) = M(ip\,ip 2 ) J — e -2*/3 2 -»/3 3 gin 2 a 

y e -l(Pl+02+03) cos Q, s j n Q, 


_ e 2j/3 2 +*ft s j n 2 a 

e 2i@ 1 +i0 3 CO g2 a 

gd/h+ZA+fo) cog a gi n a 


where the second matrix on the right-hand side of the equation above has determinant equal to 
one. □ 

As a consequence of the result here above, the determinant of M.(ipi,ip 2 ) depends only on 
the complex space spanned by ip\ and ip 2 . Therefore, in a neighborhood of a conical intersection 
between the levels A j, Aj+i we can define the following function: 

F(u) = detM(ipi, ip 2 ) (12) 


where {ipi,ip 2 } is an orthonormal basis for the sum of eigenspaces relative to the two crossing 
levels. In particular, outside the intersection we can take, for instance, i/’i = <Pj and ip 2 = 4>j+i- 
If the levels Xj,Xj + i are (locally) separated from the rest of the spectrum, the projection 
associated with the sum of the eigenspaces of the intersecting levels is continuous with respect to 
u (see Proposition A.3|, which implies that F is continuous (see |lflj). 

The following result characterizes conical intersections in terms of the conicity matrix. 


Proposition 3.2 Assume that {Ay, Ay_)_i} is a separated discrete spectrum with Aj(u) = Aj+i(u). 
Let {ipi,ip 2 } be un orthonormal basis of the eigenspace associated with the double eigenvalue. Then 
u is a conical intersection if and only if M.(ip\, ip 2 ) is nonsingular. 


Proof. Let r v (t) = u + tv, where v is a unit vector in R 3 , and let (pf,<pf +1 be the limits of 
<pj(r v (t)), <pj+i(r v (t)) as t —> 0 + (recall that the eigenfunctions <pj,(pj + \ can be chosen analytic 
along r v for t > 0). Assume that the intersection is not conical. Then for every e > 0 there is a 
unit vector v e = (vf , , r’f) such that 


d_ 

dt 


A J+ i(r Ve (t)) 




< £, 


« Ve - < e, 

i =1 

while ^ implies that b Ve := JA =1 v t (<^J e , = 0. Consider an orthogonal matrix A e having 

v e as first row. Since | det A £ | = 1, we have that 



(° 

0 

a Ve \ 

detM(</)J E ,(/>J{ 1 ) = det (A e M((/>J%^J^)) = 

det ci 

r* 

C 1 

di 

\c 2 

c 2 

d 2 J 


< C'e(o|Aj(u)| + /3) 2 , 


|2m Ve 3m(ciC2)| 
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with C, a and f3 suitable positive constants, where we have used the fact that 

for i = 1,2 and that Hi is if(u)-bounded. Thus arbitrariness of e and Lemma [l] imply that 
A4(ipi,ip 2 ) is singular for any orthonormal basis {ipi , ip 2 } of the double eigenspace. 

Let us now prove the converse statement: assume that u is a conical intersection and, by 
contradiction, that A4(ipi,ip 2 ) is singular for every orthonormal basis of the eigenspace associated 
with the double eigenvalue. We introduce the matrix 

_ (y\c((iPi,H^ 2 )) 3m ((^ 1 ,# 1 ^ 2 )) ~ (ipi,Hi'i/} 1 ))\ 

= iHc((V’i,Lf 2 '02)) 3m ((Vh,^^)) ({^ 2 ,H 2 ^ 2 ) ~(ipi,H 2 ^i)) , 

i,H 3 i> 2 )) 3m((V’i ,H 3 ijj 2 }) {(ip 2 , H 3 ip 2 ) - (i/’i, # 3 ^ 1 ))/ 

and we notice that det ip 2 ) = — 2i det i/> 2 ) so that M(4>i,4> 2 ) is singular if and only 

if A4(ipi,i/j 2 ) is. The condition of conical intersection implies that 

3 

a v = X] v i (<^+ 1 , - (<%, Hi4> J» ^ 0, Vv G R 3 . 

i =1 

Moreover, equation |8]) implies that b v = JT =1 v i{4>j 1 = 0, so that the third column of 

the matrix is never linearly dependent from the first two. In particular the matrix 

4>j+i) is singular only if the first two columns of the matrix are linearly dependent. Thus, up 
to multiplying W’W+i by a phase factor, we can always assume that £R, i = 1, 2, 3. 

Let us now fix v £ R 3 and let us call W the orthogonal complement in R 3 of the vector 
((0J,7Li0J +1 ), i?30J + i)). We have that v e W and dimW > 2. It is easy to 

prove that for every w £ W the limit basis {(j>J, ^J+i} is equal to {0J, </>J +1 }, up to exchanges 
between the two elements and up to phases. Indeed, by definition of W we have that (<j>J, [w\Hi + 
w 2 H 2 + w 3 H 3 )cj)Y +1 ) = 0, and writing 0J +1 ) T = U , (j)J +1 ) T for some U of the form ( flTj ), 
we obtain that 


0 = (0J, (w 1 H 1 + w 2 H 2 + w 3 H 3 )(/> J +1 ) 

= e~ ih {[e~ 2ipl cos 2 a - e " 2i02 sin 2 a) b w + sin a cos a a w ) , 

therefore it must be a = k^, k £ Z, since 6 W = 0 and a w 7 ^ 0. 

Let us now consider the vector 


T = «^7 +1 , H lt f,J +1 ) - (d>J, HtfJ), {<t>J +1 ,H 2 <f>J +1 ) - (<^ 7 , HrfJ), (</>7 + 1 ,ff 3 </>7 + i) - H^J)) 


which corresponds to the third column of A4((fiy,(f>Y+i)- Notice that the definition of T, up to 
a sign, does not depend on the choice of w € W, by what precedes. Since Y 3 -, the orthogonal 
complement of Y in R 3 , has dimension 2 there exists a non-zero w £ W D Y x . By definition of w 
we have that a w = w • Y = 0. 

We get a contradiction, thus the matrix M{(j) J, </>J +1 ) must be nonsingular, and therefore also 
M{(j>j , 4>j +1 ) has to. □ 

As an example, consider the Hamiltonian H( u) £ *u(3), with 

/ 0 0 0 \ / 1 i 0 \ / 0 0 i \ / -1 1 -1 \ 

H 0 = 0 0 0 , H x = -i 0 1 , H 2 = 0 1 0 , H 3 = 1 1 0 . 

\ 0 0 1 J \ 0 1—1 J \ -i 0 0 J \ -1 0 0 J 

(13) 

The Hamiltonian H(-) admits a double eigenvalue at u = 0 corresponding to the two lowest levels. 
A simple computation leads to det Al(ei, e 2 ) = —2 i where ei = (l,0,0) r ,e 2 = (0,1,0) T form a 
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basis of the double eigenspace at 0. Thus the eigenvalue intersection is conical. An example of 
conical intersection in an infinite dimensional setting is provided in Appendix |B.2| 

A peculiarity of conical intersections is that, when approaching the singularity from different 
directions, the eigenstates corresponding to the intersecting eigenvalues have different limits. The 
following proposition provides the relation between these limits. 


Proposition 3.3 Let u be a conical intersection between A j and A J+1 . Let v 0 ,v £ R 3 be two 
unit vectors, and call (p®,<p°j+i the limits as t —> 0 + of the eigenstates <pj(ro(t)), cpj + i(r 0 (t)) along 
a straight line ro(t) = u + tvo, and the limit basis along the straight line r v (t) = u + tv. 

Then, up to phases, the following relation holds: 


( <j)j \ _ f cos E e 1/3 sin E\ ( cp® \ 

VPj+ J ~ \-e il3 sinE cosE ) \f j+1 ) ’ 

where the parameters E = E(v) and /3 = /3(v) satisfy the following equations: 


tan2E(v) = (—1)^ 


(mod 27 r) 
P V = 


2|<$,.ff v $ +1 >| 

- {(t>° j+1 , H v( p 

arg(^°, H v (j>° j+1 ) + £tt, 


0 

j +1 


) 


where H v = i HiVi and £ = 0,1. 


(14) 


(15) 

(16) 


Proof. First of all, we notice that all pairs of orthonormal eigenstates of H{ u) relative to the 
degenerate eigenvalue can be obtained by the action of the group U(2) through the pair (</>°, <(>° +1 ). 
Nevertheless, we are not interested on the global phases of the states, that is we consider the 
equivalence relation (e^ipi, e*^ 2 ^) ~ (ipi, ipi) V /?i, £ R. Therefore for any unit vector v £ R 3 

we can obtain a representative of the pair (<(>J through the transformation ( |14| ), for some 

E(v) and some /3(v). This gives 


0 = (</>! ,H v q +1 ) 

= «$, H v <p° j+1 ), <$, H v( p° j+ 1 y, <$ +1 , H v <fP j+1 ) 




/ cos 2 E(v) 

J — gin 2 E(v) 
y e */3( v ) cos E(v) sin E(v) 


= (cos 2 E(v)|<$, H v( p° j+1 )\e^-^ - sin 2 E(v)|<$, H v cf ° j+1 )| e - i( ^ (v)) + 

+ cosE(v)sinE(v)((^ + 1 ,fl v ^ +1 ) - (0°, ff v </>°))), (17) 

with 7 = arg ((</>?, H v <p ° J+1 )), whenever (cp°, H v cp° j+1 ) ^ 0. 

If v = ±vo, then (cp®, H v <p® +1 ) = 0, and the conicity condition implies that E(v) = kn/2. In 
particular, /3(v) can be any real number. 

If v is not parallel to v 0 , then (cp°, H v (p°- +1 ) y 0 and /3(v) can only take the values arg (cp®, H v (p^ +1 ) + 
fir, £ = 0,1. Indeed, since the term (<pj +1 ,H v <p t j + f) — (<p®, H v (pj) is real, the imaginary part of 
the term into parenthesis of <HZf is zero: 


\{(p°,H v( p 0 j+1 }\ sin ( 7 - /3(v)) = 0. 


If we choose /3(v) = 7 , then we can prove by computation that E(v) must satisfy 


tan2E(v) = 

while if /3(v) = 7 + n, S(v) must satisfy 

tan2E(v) = — 


2|<0°,.ff v $ +1 > 


(^,F v ^)-(^ +1 ,H v ^ +1 ) 

2 |(^,U V ^ +1 )| 

(cP°H v cp<l)-(cp°H v cP°)- 


(18) 


(19) 
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Notice that the two pairs (E(v), /3(v)) and (—E(v),/?(v) + 7 r) give the same transformation in (14). 

□ 


It can be seen that not all the solutions of (15)-(16) provide the correct transformation (14), 
which, nevertheless, is easy to detect. The good solutions of (15)-(16) constitute four branches 
which are continuous with respect to v, and they can be constructed as follows. Let w(s), s £ [0, s], 
be a curve joining Vo to v such that w(s) ^ {vo, — Vo} for every s £ ( 0 , s); for conical intersections, 
it is possible to associate with such a curve a continuous solution (3(w(s)), /?(w(s))) of (15)- © 
with 3(vo) = 0 and compatible with ( fl4| . In particular, if we choose E(v) according to (18), it 
is easy to see that E(w(s)) £ [—7r/2, 0] for s £ [0, s], from which one deduces that the final value 
E(v) = E(w(s)) is independent of the chosen path and continuously depends on v. Moreover, it 
is easy to see that E(—v 0 ) = — tt/2. Similarly, one can show that /3(v) = /J(w(s)) is independent 
of the chosen path and continuous outside {v 0 , — v 0 }. Note that the fact that j3 is discontinuous 
at —Vo implies that the corresponding limit basis (</>J,^>J +1 ) has a discontinuity at — Vo, that is, 
its limit depends on the path. 

We can repeat the same argument choosing E(v) according to (19), with E(v 0 ) = 0. The other 
two continuous branches are obtained choosing the initial condition E(vq) = 7r. 


3.2 The basic control algorithm 

Let us consider the following controllability problem. 

Let H(-) satisfy (HO)-(Hl) and £(■) = {Ao(-), • ■ •, Afe(-)}. Then, given e > 0, iT, £ 
w, j £ {0, and p £ [ 0 , l ] fe+1 such that Y^i=oPi = b find T > 0 and a path 

u : [0, T\ — » w with u(0) = u s and u(T) = u-1 such that 

k 

H^( T ) _ < e, (20) 

i =o 

where iffi) is the solution of ([3]) with ip( 0) = (j>fi u s ), and r d 0 ,...,i}k G R. are some 
possibly unknown phases. 

If all levels are connected by means of conical intersection occurring at different values of the 
control, the results obtained in the previous section provide the basic elements in order to construct 
a family of control paths solving the problem here above. This can be done by taking advantage of 
the following proposition, which describes the spreading of occupation probabilities induced when 
a path in the space of controls passes through a conical intersection. 


Proposition 3.4 Let u be a conical intersection between the eigenvalues A,, A J+ i. Consider the 
curve 7 : [ 0 , 1 ] —> to defined by 


7(t) = 


u + (t - t 0 )w 0 
u T (t - r 0 )v 


t £ [ 0 , r 0 ] 
t £ [t 0 , 1 ] ’ 


for some tq £ (0,1) and some unit vectors wo,v. Then there exists C > 0 such that, for any 

£ > 0 , 

||^(l/e) - 7 r 1 e* tfl ^>j( 7 (l)) - 7 r 2 e l,? 2 0 j + i( 7 (l))|| < Csfe (21) 

where i?i, d 2 £ R, if(-) is the solution of equation § with = (pj( 7 ( 0 )) corresponding to the 
control u : [ 0 , 1 /e] —l w defined by u(f) = 7 (et), 


7 Ti = | cos (E(v)) j, 7 T 2 = |sin(E(v))|, 


and E(-) is 
H(-w 0 ) = 
respectively. 


the only solution of equation PL such that E(v) £ (— 7 t/ 2 , 0 ] for v 7 ^ -wo, and 
— 7 t/ 2 , where the limit basis in (18) is given by the limits </>j( 7 ('r c /)),</>j+i( 7 (T c j")), 
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Proof. We consider the Hamiltonian H(u(t)), t £ [0,1/e]. Since the control function u (•) is 
not C 2 at the singularity, we cannot directly apply the adiabatic theorem. Instead, we consider 
separately the evolution on the two subintervals (in time t) [0,ro/e] and [lo/e, 1/e]. 

Since the eigenstates </j(u(f)),/ 7 + i(u(f)) are piecewise C 1 , we can apply [30] Corollary 2.5] 
and obtain that there exists a phase 6\ (depending on e) such that 

||^(r 0 /e) - e l 01 0 j( 7 (r o _ ))|| < C\[e, 

for some constant C > 0. By Proposition |3.3[ this implies that 

HV’^o/e) - e lf>1 (cosH(v) </> J ( 7 (r 0 +)) - e" l/3(v) sinH(v) ^•+i( 7 ('r 0 + ))) II < c Ve, 

with S(v) as in the statement of the proposition and /3(v) = arg(</°, H v (/>^ +1 ). 

Applying PHI Corollary 2.5] also in the time interval (ro/e, 1/e], we conclude that there exists 
two phases di and $ 2 (depending on e) such that 

||^(l/e) - (e 101 cosH(v) <^( 7 ( 1 )) - e Ilt>2 sinH(v) ^+ 1 ( 7 ( 1 ))) II < C'y/e, 

for some constant C > 0. □ 

For control purposes, it is interesting to consider the case in which the initial probability is 
concentrated in the first level, the final occupation probabilities p\ and p 2 are prescribed, and we 
want to determine a path that induces the desired transition. For a given line reaching the conical 
intersection, the outward directions that provide the required spreading of probability are given 
in the following proposition. The proof follows from simple computations and is thus omitted. 


Proposition 3.5 Let u be a conical intersection between the eigenvalues and let 717 ,7r 2 

be positive constants such that nf + n 2 = 1. Consider the line r(t) = u + (t — to)wo, t £ [0, To], 
for a unit vector wo £ M 3 and some tq £ (0,1), and set (fi® = ^>j( , y{ T o)) and = ‘/’j+iW'hT))- 
Then the locus formed by the directions v £ R 3 that give rise to transformation ( 14) with 


7Ti = | cos^(v)| and 7 r 2 = | sins(v)| is given by the following expression whenever ^ {0,1} 

where K = { (x, y, z) £ I 3 : \Jx 1 + y 2 = Cz} and 

7T17T2 


( 22 ) 


C = 


^ ■ 


Otherwise, 1/717 =0 then v = Wq and if 717 = 1 then v = — Wq . 


The controllability problem presented at the beginning of this section can be solved taking advan¬ 
tage of the results shown above. The strategy consists in constructing a piecewise C 2 path joining 
u s with tF that passes through the conical intersections ft, between the j- th and the (j + l)-th lev¬ 
els, j = 0,..., k — 1, and avoids any other degeneracy point. The tangent directions at the conical 
intersection are chosen according to the probability weights p 2 , as explained in Proposition 

When we are far from all the conical intersections, we approximate the evolution with that of 
the adiabatic Hamiltonian 


3.5 


h a (r) = H{i{t)) - «£^P/(r)Pi(r) - ieP- L (r)P- L (r), 


(23) 


1=0 


where P;(t) is the spectral projection onto the eigenspace relative to Aj(y(r)) and P _l (t) = 


id — 2Z)z = 0 Pi{ T )- The evolution associated with (23) conserves the occupation probabilities relative 
to each energy level A;, l = 0,..., k, then so does the evolution of H(p/{t)), with an approximation 
of order e (see Remark [I]). 


11 







In a neighborhood of the conical intersection between Xj and A 7 + 1 , j = 0, ...,k — 1, we 
decouple the evolution inside the band of the two intersecting levels from that relative to the rest 
of the spectrum, that is we approximate the dynamics with the ones associated with the adiabatic 
Hamiltonian 


k 

K(t) = - ieP(T)P(T) - is ^ Pi{t)Pi{t) - *eP J -(r)P- L (r) (24) 

1=0 

1 ^ 3 , 3 +1 


where P(r) is the spectral projection relative to {Aj( 7 (r)), Aj+iWi"))}- The evolution associated 
with (24) conserves the occupation probability relative to the band {Ay, Aj_|_i}, the ones relative 
to all other energy levels in {Aq, ..., A*,} and the one associated with the rest of the spectrum. 


The evolution given by (24) inside the band {Aj, A J+ i} is described in [301 Corollary 2.5] and in 
Proposition |3.4| 

Let us now explicitly determine the path in a neighborhood of each conical intersection. With¬ 
out loss of generality, we assume that ^( 0 ) = <^o(u s ); in the other cases a path can be obtained 
similarly. We construct a (piecewise C 2 ) path 7 : [0,1] —)• w as follows. 

We set 7 ( 0 ) = 11 s and, for some 0 < r 0 < 1 we choose 7 |[o,r 0 ] in such a way that 7 (r 0 ) = Uo 
and all the eigenvalues A;( 7 (r)) are simple for every l = 0,..., k and r £ [0, r 0 ). Moreover, 7(3 
is chosen tangent to a segment in a neighborhood of Uo, and we call Vjj" the tangent direction 
of j(-) at t = Tq . In a neighborhood of u 0 and for r > r 0 , 7 ( 3 , is cho sen to be tangent to v[j" 


'o 

where the outward direction v, 
tt2 = y/P 1 


Pi 


+ ; 
0 


is chosen according to Proposition 3.5 with Wq = v 0 , 717 = po, 


The rest of the path is constructed recursively. Assume that the path has been defined up to 
time Tj- 1 , for some j = 1 ,..., k — 1 , with 7 ( 77 - 1 ) = Uj_i, and that an outward direction 
has been selected. Then choose 77 £ ( 77 -- 1 ,1) and the smooth path 7 |( T ._i, T .](-) joining Uj-i with 
Uj such that all the eigenvalues A;( 7 (t)), are simple for every l = 0 ,..., k and for r G (tj_ 1 , 77 ), 
in a neighborhood of Uj_! the curve is tangent to v+ and in a neighborhood of Uj the curve is 


tangent to some segment that we call 


apply Proposition 


3.5 


with wq = v 


3 ' 

7Ti = 


To choose the outward direction of y(-) at u 7 
-HII- and 7 T 2 = y/l — Trf. 


we 


The last arc defined on ( t^-i, 1] is simply constructed by joining u^._i with u-1, taking care of 
choosing y(-) tangent to the outward direction in a neighborhood of Ufc_i. 

To avoid highly non-homogeneous parameterizations, the path 7(3 can be reparameterized by 
arc-length. 

Let us now reparameterize the time setting t = r/e, for some small positive e. The adiabatic 
theorem and Proposition |3.4| lead to the estimate 


||V’( 1 /e) - < Cyfe, 


1=0 


for some do,..., G R and some C > 0 depending on the path y(-) and on the gaps in the 
spectrum. The geometric construction of the path j(-) is represented in Figure [ 2 ] 


4 An improvement of the efficiency of the algorithm: the 
non-mixing field 

The problem of reducing transition times is quite important in quantum control, in particular for 
the need of reducing decoherence. Therefore, in this section we will show how to do it by means 
of some special paths in the space of controls that eliminate, in the adiabatic approximations, 
the error coming from the intersection of the eigenvalues, so that the only error induced by the 
application of the adiabatic theorem comes from the gap among the remaining eigenvalues and 
those among these eigenvalues and the band of the intersecting ones. 
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Figure 2: Construction of the path 7 (-). The corners at the points U;, corresponding to conical 
intersections, are chosen in such a way that they induce the desired spreading. 


Let us consider a pair {A j, of eigenvalues, which are simple and separated from the rest 

of the spectrum, according to assumption (HI), in a certain open set ui C R 3 , except for a point 
u £ u, where they have a conical intersection. We are interested in the dynamics inside the 
subspace P u ?l, where P u denotes the projection associated with the two levels {Aj(u), A J+ i(u)} 
for u € to: we know that, under adiabatic evolution, the dynamics are described by the effective 
Hamiltonian (101. To improve the precision of the result, the idea is to cancel the off-diagonal 
terms in the effective Hamiltonian, which are responsible of the error of order yfe the estimates 
given in EDI Corollary 2.5] and in Proposition |3.5| In order to do that, we choose some special 
trajectories in ui along which the term (fj . is null, that is, thanks to (J7|, we look for curves 

y(-) in the space of controls that satisfy the equation 


(7i(t)-ffi +7 2 (i)#2 +'y 3 (t)H 3 )cl> j+1 (t)) = 0. (25) 

We denote the first column of the conicity matrix AI( - i/’i) V’2) by 

"iWb, 1P2) = ((Vh, #1^2), {tpi, H 2 %f 2 ), (1P1, H 3 ip 2 )) T , 


and its components (^i, Hiip 2 ) as m;. It follows by definition that the real vector 

Xyifi,^) =- Yi - 

= (3m(m 2 m 3 ),3m(m3m'l),3m(mim 2 )) T , 

where x denotes the cross product, is orthogonal to both m(ifi,if 2 ) and m*(if\,if 2 )- 


(26) 


Remark 2 Let us remark that the vector X(ipi, ip 2 ) is invariant under phase changes in the 
argument, that is X(ipi,ip 2 ) = X(e l ^ 1 ifi,e l ^ 2 if 2 ). Notice however that Xfifi,ifi 2 ) = —X(ip 2 ,ipi). 


Definition 4.1 The vector field 


X P (u) = X(<t>j(u),<l> j+1 (u)), 


(27) 


defined in w \ {u}, is called the non-mixing field associated with the conical intersection u. 

The non-mixing field is smooth in its domain of definition. From 0 and ( p6| ), we have (ct>j , <fij+i) = 
0 along its integral curves. Moreover, a simple computation leads to 


/{if 2 ,H 1 ijj 2 ) - (4> 1 ,H 1 'ip 1 )\ 1 

Xfif 1 ,^ 2 ) ■ {i>2,H 2 -ip 2 } ~ (%/j 1 ,H 2 ip 1 ) = — detM{ipi,ip2), ( 28 ) 

\(^2,H 3 ip 2 ) - (i/’i^V’i)/ 1 


which implies the following result. 
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Proposition 4.2 Let y(-) C u>\ {u} be an integral curve of the non-mixing field. Then 

- Aj(7 (t))] = ^-.Ffiy(t)). 

In particular, all the integral curves of the non-mixing field starting from a punctured neighborhood 
of the conical intersection reach it in finite time (up to a time reversal). 

Without loss of generality, we can assume that u = 0. Denote every u £ I 3 by u = pv, 
where p is its module and v the versor. Let us call p(p,v) = m((j)j(p, v), <j>j + i(p, v)), where 
<f>j(p, v),<j)j + i(p, v) denotes a choice of the eigenstates relative respectively to A j(pv) and Aj+i(pv) 
(extended by continuity for p = 0). We remark that p(p, v) is defined up to a phase, therefore its 
module and the quantity 

p(p, v) x/i'(/),v) = 2 i 3m(p(p, v)) x e(p(p, v)) (29) 

are well defined. Finally, set X^(p, v) = X P {pw). Notice that, for different values of v, X^(p,w) 
has a priori different limits as p —> 0. 

The following estimates hold. 

Lemma 2 Assume that u = 0 is a conical intersection. Then the inequality 

KP> V ) • v| < Cp (30) 

holds in a neighborhood of 0 for some constant C > 0 uniform with respect to v. 

Proof. Without loss of generality, we assume that the double eigenvalue is equal to 0. Then 
pp(p,v) • v = p(<j>j(p,v), {viHx +v 2 H 2 + v 3 H 3 )<f> j+ 1 (p,v)) 

= -(<Aj(A v ),^o0j+i(p,v)) = -(<Mp,v) - Po<(>j(p,v),-ffo0m(p- v ))- 

Assumption (HO) implies that 


3 

\\Ho(f>j+i{p, v)|| = \\\j + i(pv)(f)j + i(p, v) -^2pViHi^ j+1 (p,v)\\ 

i =1 
3 

< |Aj+i(p,v)| +P^ ||^</>j+i(p,v)|| 

i-1 

< cp, 

for some c > 0, locally around the intersection. By smoothness of the projection, we get that 

|(<MP,v) - P 0 <l>j{p,v),Ho<f> j+1 (p,v))\ < Cp 2 , 

for a suitable C' > 0, then we get the thesis. □ 

We are now ready to prove the following result, which provides some information on the 
behavior of the trajectories of the non-mixing field. 

Proposition 4.3 With the notations introduced above and for p small enough, there exist three 
constants C\,c 2 ,C 3 > 0 such that c\ < \p\ < c 2 and |v| < C 3 along the trajectories of the non-mixing 
field. 

Proof. Direct computations lead to the equations 

P = X,i ■ v, 

v = - {X^ - (X„ • v)v). 
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The upper bound for |p| comes easily from the -ff (u)-boundedness of Hi for every i. 

From ( [28| and p 9 | , we get that \X^\ = |3m(p(p,v)) x 91e(/i(p, v))| > c for some c > 

0 in a neighborhood of the singularity. An immediate consequence is that the three vectors 
91e(p(p, v)), 3m(p(p, v)) and X ^(p, v) are three linearly independent vectors in R 3 . Then we define 
the corresponding three unit vectors Zi(p, v) = 9te(/i(p, v))/|91e(p(p, v))|, z 2 (p, v) = 3m(p(p, v))/|3m(p(p, v))| 
and w(p, v) = X^/\X^\. We can then decompose v in this basis as 

v = a(p,v) Zl (p,v) + 6 (p,v)z 2 (p,v) + 7 (p,v)w(p,v). 


Notice that there exists some S < 1 such that |zi(p, v) • z 2 (p, v)| < S for every (p, v). Estimating 
the scalar product of v with zi and z 2 , using ( |30| ), we obtain that 

Kp,v) + 6 (p,v)z!(p,v) • z 2 (p,v)| < Cp | b(p,v) + a(p, v) z^p, v) • z 2 (p,v)| < Cp, 


that easily leads to |a(p, v)| < Cp and \b(p, v)| < Cp , for some uniform constant C. 
From (30) and the fact that |v| = 1 we get that 


|w(p, v) • v| > | 7 (p, v)w(p, v) • v| 

= l(v- (a(p,v)zi(p,v) + &(p,v)z 2 (p,v))) -v| 

= |1 - |a(p,v)zi(p,v) + &(p,v)z 2 (p,v)| 2 | > 1 - 4Cp 2 . 


This implies that 

|w(p, v) - (w(p, v) • v)v | 2 = |v - (w(p, v) • v)w(p, v )| 2 = 1 - (v • w(p, v )) 2 < 8 Cp 2 . 
The thesis comes from the fact that 


\{Xp - {Xu, • v) v)| = |A M ||w(p,v) - (w(p,v) • v)v| and | Xy, ■ v| = \X^\ |w(p, v) 


□ 

The following proposition is a generalization of m Proposition 5.9] in the three dimensional 
case. The proof follows the same lines, thanks to Proposition |4.3[ and is thus omitted. 

Proposition 4.4 For every unit vector v in 
Xp with 7(0) = 0, p > 0, such that 


there exists an integral curve 7 : [— 77 , 0 ] —> co of 


lim 

t —>-0“ 


7 (*) 

\m\\ 


Thanks to Proposition the integral curves of the non-mixing field are C 1 up to the singular¬ 


ity included. In particular, they satisfy the hypotheses of Proposition A.4 with k = 1, so that the 
projections Pj{ u) and P, + i(u) on the eigenspaces relative to the intersecting eigenvalues are C 1 
along the integral curves of the non-mixing field outside the singularity, and can be continuously 
extended at the singularity. On the other hand, P u is C 1 along such curves, singularity included. 

We remark moreover that, if 7 : [to,ti] —> R 3 is an integral curve of the non-mixing field such 
that Aj( 7 (t)) 7 ^ Aj+i( 7 (f)) for t £ [to,ti), by definition of the non-mixing held, it holds 


Pj( 7 (t))P j+1 ( 7 (t)) = 0 P j+1 ( 7 (t))Pj( 7 (t)) = 0 


V t £ [fo,fi). 


(31) 


We have the following. 

Proposition 4.5 Along every integral curve of the non-mixing field, there is a choice of the 
eigenstates relative to the intersecting eigenvalues which is C 1 up to the singularity included. 

Proof. Let I = [—T, 0], and let 7 : I —> R 3 be an integral curve of Xp such that 7 ( 0 ) = u is a 
conical intersection between A j and A ;+1 . Outside the singularity, the eigenstates are well defined, 
up to a phase. To hx the phase, we set 

„ (t) = mmi , 

IIPiMWII 
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where if is an eigenstate of H( u) relative to Aj(u) = Aj+i(u) such that lim t _>.o Pji'lit))^ ^ 0- 
Possibly reducing T, we can assume that ^ 0 on the whole I, thus ipj(t) is a normalized 

eigenstate of H(j(t)) relative to Xj^ft)). Since ||P/( 7 (i))V’|| ^ 0 , in order to prove that ipj{t) is 
C 1 it is enough to prove that Pj^ft))^ is. 


Since P,( 7 (t)) + P J+ i( 7 (t)) = P 7 p) for t G [-T,0), and by ([31]), we get that Pj{^{t))Pj{l{t)) = 
Pi( 7 (i))P 7 ( t ), an d therefore 




WPMWW 


which has limit for t —> 0 


Let us notice that P 1 ^Pj(^{t))ip = Pj('y(t))ip. This, together with (311, implies that 


where the right hand side has limit for t —> 0 _ . 

We can repeat the same procedure to show that there is a choice for ipj+i (f) such that ipj + ±(t) 
has limit for t —¥ 0 _ . □ 

As an immediate consequence of the above proposition, we get the following corollary. 

Corollary 4.6 In a neighborhood of a conical intersection, the integral curves of its associated 
non-mixing field are C 2 up to the singularity included. 

The regularity results proved earlier can be improved as shown below. 

Proposition 4.7 In a neighborhood of a conical intersection, the integral curves of its associated 
non-mixing field are C°° up to the singularity included. In particular, we can choose C°° eigenstates 
ifj, ifj+ 1 along such a curve, up to the singularity included. 

Proof. Let 7 : [— T, 0] —» R 3 be an integral curve of Xp such that 7(0) = u and |7(t)| > 0 
for every t G [—T, 0] (this is true up to choosing T sufficiently small), and define the eigenstates 
1 (-) as in Proposition 4.5 We prove by induction on n that, if the curve is C n {[—T, 0]), 
then also ipj(') and ipj+ l(-) are C n ([—T, 0]). In particular, this last fact implies that the integral 
curves of the non-mixing field are C n+1 ([—T, 0]). 

In Proposition |4.5| the claim is proved for n = 1. Assume that it is true up to n — 1, with 
n > 1, that is ipj(-) and i/jj+ 1 (-) are C n ~ 1 ([—T, 0]), and then the integral curves of the non-mixing 
held are C n ([—T, 0]). In particular, we know that P 7 (.) G C”([—T, 0]), and Pj(7(-)), Pj+i(7(-)) G 
C n -\[-T, 0]). 

Differentiating n times the identity = P 1 (p j ifj(t), we obtain that for every t G [— T, 0) 

ip^ (t) is a linear combination of the terms ^ (' t ), for l = 0,..., n. By inductive hypoth¬ 

esis, all the terms relative to l > 1 are known to be continuous on [—T, 0]. Then we have to prove 
the continuity of P 7 ( t )'i/’]”' l (t) at t = 0, which is equivalent to the continuity of P 7 ( t )Pj"' l ( 7 (t))^>. 
We can develop P^ t )Pj n \'y(t)) as 


ppffi) = p.p)") + P J+ 1 P^ n) 


= P,-P<"> - P,Pt\ + P, 


J'j + l -T*'j+lP; 


( n ) 


where we omitted to precise that all projections are evaluated along 7 for simplicity of nota¬ 
tion. Notice that the term PjP ' n ^ is continuous on the closed interval. By (31), it follows 

d t (Pj(. r y(t))Pj+i('y(t))') = M G [—T, 0), then we can write Pj( 7 (t))Pj+i( 7 (t)) as a linear 


dt n ~ 
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Figure 3: Trajectories of the non-mixing field for the Hamiltonian corresponding to (13) 


combination of the terms Pj l \'y(t))Pj™f l \'y{t)), for l = 1,..., n — 1. Since all these terms are 
continuous on [— T, 0] by induction hypothesis, it follows that also Pj('y(t))Pj T ^ 1 ('y(t)) is. 

Analogous computations prove the same for Pj + i("/(t))Pj n \"/{t)), therefore also PP- n ' > can be 

defined continuously on [—T, 0]. In particular, P^^rp^ft) is continuous on [—T, 0], and we get 
the thesis. The same argument proves the smoothness of ipj+i(t). □ 


Figure [3] describes the behavior of the trajectories of the non-mixing field relative to the first 
two eigenvalues of the Hamiltonian H{ u) £ iu(3) defined by (13). Consistently with the results 
shown above, the flow corresponding to the non-mixing vector field allows to identify two conical 
intersections, one of them being the origin, among the two levels. In particular, the trajectories 
converge or diverge from them, locally. 

To conclude this section, we present below a result providing some information on the structural 
stability of conical intersections based on the properties of the non-mixing fields. 


Theorem 4.8 Assume that H{ u) = Hq + UiHi + U2H2 + U3H3 satisfies (HO)-(HI) and let u 
be a conical intersection for H (u) between the eigenvalues A j and Aj+i belonging to the separated- 
discrete spectrum S. Then for every e > 0 there exists S > 0 such that, if H( u) = Hq + u\H\ + 
U2H2 + U3H3 satisfies (HO)-(Hl) and 


3 

Ei 

i=0 


Hi - H 


H 0 


<5, 


(32) 


then the operator H( u) admits a conical intersection of eigenvalues at u, with |u — u| < e. 

Proof. First of all, by equivalence of all norms || • ||#( u ), without loss of generality we can 
assume that u = 0. We notice that our assumptions guarantee that in a neighborhood of the 
conical intersection the eigenvalues A j and A J+ i are well separated from the rest of the spectrum. 
Continuous dependence of the eigenvalues with respect to perturbations of the Hamiltonian (see 
Lemma [ 5 ]) ensures that, if 6 is small, then H{-) admits two eigenvalues A,-, A-/ + -| close to A ? , A J+ i. 
Moreover {A,, A )+ i} is separated from the rest of the spectrum, locally around u. 

From the conicity of the intersection between A j and A J+1 , there exists e > 0 small enough 
such that |F(u)| > c for some c > 0 on B( u, e) and moreover, by Proposition 4.3 the vector 
field Xp (up to a global sign) points inside the ball £?(u, e) at every point of its boundary. If 
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S is small enough then Xj ^ A J+ i on dB( u, e) and the gap between the two eigenvalues can be 
assumed to be of order e. Therefore we can define the conicity matrix Ad associated with H(-) and 
the function F( u) = det Ad(^(u),^- + i(u)), where {{^(u),{7+1(11)} is an orthonormal basis for 
the sum of eigenspaces relative to {Aj(u), Aj+i(u)}. Since the conicity matrix varies continuously 
with respect to the control operators and the vectors on which it is evaluated, we can take 6 small 
enough such that F(u) > c/2 on B(u, e). This allows us to define, whenever Xj ^ Aj+i, the 
non-mixing field Xp associated with H(-) and corresponding to the band {Ay, Ay_(_i}; thanks to 
up to a time reversal the time derivative of A^+i — Xj along the integral curves of 

if S is small enough, then Xp points inside B{ u, e) at 


Proposition 


4.2 


A.2 


Xp is smaller than —c/4. By Corollary 
every point of dB(u, e). 

Any trajectory {(-j of Xp starting from B(u,s) remains inside B(u,e) in its interval of defi¬ 
nition and reaches in finite time a point u corresponding to a double eigenvalue Aj(u) = Aj+i(u). 
The conclusion follows from Proposition |3.2| □ 


4.1 A spread controllability result 

In this section we show how the curves tangent to the non-mixing field allow to improve the 
performances of the control algorithm presented in Section [372] The result is stated here below. 


Theorem 4.9 Let H( u) = Ho+uiHi+U 2 H 2 +U 3 H 3 satisfy hypotheses (HO)-(Hl). Assume that 
there exist conical intersections Uj G cu, j = 0,..., k — 1, between the eigenvalues Xj,Xj+i, with 
A;(uj) simple if l ^ j, j + 1. Then, for every u s and cd such that the eigenvalues A;, l = 0,..., k 
are non degenerate at u s and re , for every </> G {</>o(u s ),..., </>fc(u s )}, and p G [0, l ] fc+1 such that 


J2j=oPj = b there exist C > 0 and a continuous control y(-) : [0,1] —> 

7 ( 1 ) = ib, such that for every e > 0 

k 

HV’li/e) - II b Ce ’ 

j=0 

where tp(-) is the solution of Q) with ip( 0) = <f>, u(f) = 7 (et), and do, ■ ■ ■ ,dk G 
depending on e and 7 . 


with 7 ( 0 ) = u s and 


(33) 


are some phases 


Proof. The strategy is analogous to the one presented in Section [312] and consists in constructing 
a piecewise smooth path joining u s with cd that passes through all the conical intersection; in 
particular, we assume that the path 7 : [0,1] —)• w satisfies 7(0) = u s , 7(1) = rd and 7(77) = 
Uj, j = 0, ...,k — 1, for some 0 < r 0 < • ■ ■ < T 7 _ 1 < 1. The only difference concerns the 
choice of the paths in the neighborhoods of each conical intersection: indeed, in these regions our 
paths are chosen to be tangent to the non-mixing field. At the intersection, the inner and outer 
directions are selected according to Proposition |3.5[ as explained in Section [3~2| and the existence 


of corresponding trajectories tangent to the non-mixing held is guaranteed by Proposition 4.4 
Let us show that, given such a curve 7 (-), the adiabatic approximations lead to the esti¬ 


mate (331. 


Far from all the conical intersections, we approximate the evolution with that of the adiabatic 


Hamiltonian (231, which conserves the occupation probabilities relative to each energy level in the 


separated discrete spectrum. 

On the other hand, in a neighborhood of the conical intersection between Xj and Aj+i, 
j = 0,..., k — 1, we approximate the dynamics with the ones associated with the adiabatic Hamil¬ 


tonian (24). 


We show the details concerning the passage through Uo, and the others can be treated anal¬ 
ogously. Since the path is tangent to the non-mixing field, we can apply Theorem 2.4 in order 
to study the evolution inside the space P u (t)'H. For r in a left neighborhood of r 0 , we can then 
construct the effective Hamiltonian and its associated evolution operator which is diagonal, 
which implies that there exists a phase 0q (depending on e) such that 

||V>(to/£) - e l6 °(foixijo ))|| < C 0 £, 
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where ip(-) is the solution of equation ([3]) with ^(0) = </o(7(0)) corresponding to the control 
u(i) = 7 (et), defined on [ 0 , 1 /e]. 

By Proposition |3.3[ this implies that 

II if>(ro/e) - e lf>0 (cosE(v) ~ e~ t/S(v) sinH(v) || < C 0 e, 


with S(v) and /3(v) satisfying equations (15)-(16), and v is the outer direction. 

Since the effective Hamiltonian is diagonal also for r belonging to a right neighborhood of To, 
we conclude that there exist two phases ao and u\ (depending on r and e) such that 


Wipir/e) - e ia °p 0 Mli T )) - <M7(t))II < C 0 e. 

Since analogous estimates hold on each passage through a conical intersection and outside the 
corresponding neighborhoods the theorem is proved. □ 


5 Final remarks 


The control strategy presented in the last section (Theorem |4.9[ ) highlights the role played by the 
integral curves of the non-mixing field to obtain controllability results with an approximation of 
order e on time intervals of order 1/e even in neighborhoods of conical intersections, while the 
classical theory would guarantee an error of order yfe. 

As proved for the two-inputs case in m . also for the three-inputs case it is possible to see that 
the 1-jet of the control path at the conical intersection determines the target probabilities, while 
the 2-jet is responsible of the error introduced by the adiabatic approximation. In other words, C 2 
approximations (in the sense precised just below) of the integral curves of the non-mixing field also 
ensure an adiabatic approximation of order e. Indeed, let us consider an arc-length parametrized 
curve 7 p (r), r small, tangent to the non-mixing field, that reaches a conical intersection between 
A j and Xj+\ at time r = 0, and let 7 be a C 3 curve such that |7(r) — 7 P (t)| < Ct 3 , for r small 
enough and some positive constant C. In particular, Proposition |A.4| implies that there exists 
a constant C' > 0 such that |</z( 7 (t)) — <M7 P (t))| < C't 2 and |</>z( 7 (t)) — </>z( 7 P (t))| < C't, 
l = j, j + 1, where as usual <j> z(u) denotes the eigenstate relative to A/(u) S E(u), evaluated at u. 
Let t = t/e and consider the effective Hamiltonians evaluated along the two curves. It is then easy 
to see by simple computations that there is a constant C" > 0 such that the difference between 
the two effective Hamiltonians is less or equal than C”e 2 t. This term, integrated along a time 
interval of order 1/e, gives a difference of order e. 


An interesting controllability problem alternative to the one introduced in Section 3.2 aims 
at sending (approximately) an initial state ip s = ]Cj=o c i < A?( uS ) 1° a final state concentrated in a 
single energy level. This problem may appear completely equivalent to the previous one, but is 
actually more delicate, due to the presence of relative phases among the levels in the initial state. 
Indeed, a natural way to induce the desired transition would be to run backward in time one of 
the paths that produces any state with the same probability distribution as il> s starting from the 
concentrated state, constructed as in Section [3.2| and Theorem |4.9| However, a simple computation 
shows that, at each passage through a conical intersection, the components corresponding to the 
intersecting eigenvalues recombine in a concentrated state only if their relative phase coincides 
with the one induced by the unitary transformation of the limit basis (that is, the phase /3(v) of 
Proposition 3.3). 

There are several strategies that in principle could overcome this issue. For instance, when 
following the given path backward in time and before reaching a conical intersection, it is always 
possible to stop for a certain time period at some point in the control space in order to control 
the relative phase between the two intersecting levels. Note that in the adiabatic evolution the 
effective stopping periods depend on the chosen speed e, while the geometric path in the space of 
controls does not depend on it. An alternative strategy consists in exploiting the non-uniqueness 

it is indeed possible to see that, for a path reaching a 


property underlined in Proposition 3.5 
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conical intersection and for any superposition of the two intersecting levels, there always exists a 
choice of outward direction allowing to concentrate the probability on a single level. Nevertheless, 
since the path is determined taking into account the dynamical phases, its construction depends 
on the time parameterization. 

It is clear that the main drawback of the methods here above is that they rely on the com¬ 
putation of dynamical phases, which comes from the integration of the energy on intervals whose 
length is of order 1/e and which are very sensitive to changes in the speed e. This compromises 
the constructiveness of the algorithm. 

Let us finally mention that, to further improve the above controllability property, under as¬ 
sumptions (HO)-(Hl) it is possible to modify the strategy, again with non-constructive arguments 
(for instance, by exploiting the rational independence of the gaps between the eigenvalues under¬ 
lined in |12[ Lemma 14]) in order to approximate at the final point not only any given choice of 
probability weights associated with E, but also any choice of the corresponding phases. 

Summing up, the constructions above combined with the control algorithm in Section |3.2| 
provide the following approximate controllability result: 

under assumptions (HO)-(Hl), assuming that all energy levels in the separated dis¬ 
crete spectrum are connected through conical intersections, and for any given initial 
and target states 0 s , 0^ distributed in E and e > 0, there exists a control input steering 
the system from 0 s to a final state whose distance from ij}f is less than e. 

It is opinion of the authors that all the results here above still hold in a non-linear sufficiently 
smooth setting, that is for Hamiltonians of the form H( u) whose derivatives with respect to the 
parameter u are IL(uo)-small up to a suitable order and under hypothesis (HI). This case is 
interesting, since it covers relevant physical models, such as those described by Hamiltonians with 
controlled electromagnetic potentials. This topic will be the subject of further studies by the 
authors. 


A Regularity properties 

Let TL be a complex separable Hilbert space; all operators in the following are assumed to be 
operators on TL. In this section we derive some regularity results on the eigenvalues and the 
eigenstates of self-adjoint operators with respect to the norm defined in |2|. The results here 
below - partially already known in literature (see for instance m) - are proved by classical means. 

In the following, p(A) denotes the resolvent set of the operator A and cr(A) its spectrum. The 
resolvent of A in 0 e p(A) is denoted by R(A,<j) = (A — 0id) -1 ; we recall that it is a bounded 
linear operator that maps TL into D(A), and that, given two self-adjoint operators A±, A 2 with the 
same domain, their resolvents satisfy the Second Resolvent Identity 

R(A 2 , C) - R(A U C) = R(Au C)(Ai - A 2 )R(A 2 ,(). (34) 

First of all let us state the following technical lemma, which will be largely used in the following. 
Its proof easily comes from the definition of || • ||a and is thus omitted. 

Lemma 3 Let A, B be self-adjoint operators with B A-bounded and 0 € p(A). Then the following 
inequality holds: 

II BR{A, Oil < (l + (ICI + 1) II R(A, Oil) ||B|U- (35) 

The following result shows that the resolvent set for a self-adjoint operator A enjoys some 
continuity properties with respect to small perturbation in the space C(T>(A),TL). The proof 
follows from the definition of resolvent set p{A) and properties of the resolvent R(A, Q. 

Lemma 4 Let A\ be a self-adjoint operator, and assume [Ci , C 2 ] C p(A\) for some real 0 < O- 
Then there exists a 6 > 0 such that if ||Ai — R 2 IU 1 < <5, then A\ and A 2 have the same domain 
and [Ci, C 2 ] C p{A 2 ). Moreover, the inequality 

11^2,0-^,011 <SC (36) 
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holds on [£ 1 ,( 2 ] for some constant C depending on C1X2 and A±. 


Proof. Let ( € p{A\). If the operator id + ( A 2 — Ai)i?(Ai,£) is invertible then the resolvent 
R(A 2 ,() is we ll defined and bounded, and satisfies 

R(A 2 ,0 = R(A 1 ,o(id + (A 2 - AJRiAuOy 1 . 


Thus the thesis follows once proved that, for every ( £ [Ci, C 2 ] , we have ||(v4 2 — Ai)R(Ai, £)|| < SC' 
for ||Ai — ^4_2II^ 4.1 < S with S small enough and for some C' > 0. This fact is a consequence of the 
uniform boundedness of ||i?(Ai, £)|| on [Ci j C 2 ] (see TS]) and from (35), and the thesis holds with 

C = 2 max Ce[fl , C2] ||fl(A lt C)|| (l + (ICI + 1) P(^i, OlO • □ 


Let A £ cr(A) be an eigenvalue of the self-adjoint operator A. For every positively-oriented 
closed path T C C encircling A, and not encircling any other element in cr(A), the projection P 
onto the eigenspace relative to A is given by 


P= -(2-Ki)- 1 £ R(A, C) d(. 


Proposition A.l Let A± be a self-adjoint operator, and let X be a simple eigenvalue of A\ such 
that er(Ai) D [A — g, A + g] = {A} for some g > 0. Then for every e > 0 there exists a 5 > 0 
depending on g and on |A| such that if \\Ai — ^ 2 1|^4i < S, then 

i) 17 (^ 2 ) n [A — g, A + g] is made of only one point p, which is a simple eigenvalue of A 2 ; 

ii) Calling P^ 1 the projection onto the eigenspace of A\ relative to A and P A2 the projection 
onto the eigenspace of A 2 relative to p, it holds 


I P^-Pu 2 


< e. 


Proof. From preceding lemma, for every S > 0 small enough, if \\Ai — -A 2 1|^4i < S, then a(A 2 ) D 
[A — g, A + g] is contained in the interval (A — g, A + g). Let T be the circle in the complex plane 
of radius g centered at A, and consider the projection P A2 = — (27n) _1 f r R(A 2 ,() d(. From p6| 
we obtain that 


tjA\ 
r X 



< (27T)- 1 11^,0-^2,011 d( 

< gCS 


that is, || P Al — P A2 \\ < 1 for i5 small enough, which easily implies that dimRange(P' 42 ) = 
dim Range(P/ 1 ) = 1, therefore a(A 2 ) n [A — < 7 , A + g\ contains exactly one spectral point p, which 
is a simple eigenvalue for A 2 (see [251 Theorem XII. 6 ]). 

In particular, P A2 = P A2 is the projection on the eigenspace relative to p , and satisfies ii) for 
S small enough. □ 


Corollary 

that if ||.Ai 


A.2 Under the hypothesis of Proposition \Afl\ for every e > 0 there exists a 8 > 0 such 
- A 2 IIyq < 5, then 


Hi 1 -<^ 2 II < 


where <p Al and (f> A2 denote respectively the eigenstate of A\ corresponding to A and the eigenstate 
of A 2 corresponding to p (normalized and with a particular choice for the global phases). 


The next result provides an estimate concerning regularity properties of the eigenvalues. 
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Lemma 5 Let A\ be a self-adjoint operator such that er(Ai) n I is discrete and without finite 
accumulation points for some open, possibly unbounded, interval I. If 6 > 0 is small enough and 
A 2 is a self-adjoint operator satisfying ||A 2 —AiH^ < S, then the eigenvalues of A 2 contained in I 
are close to those of Ai, in the following sense. Up to appropriately indexing on a subset of'L the 
eigenvalues (counted with multiplicity) in cr(Aj) Dl, for j = 1,2, and denoting them with pi(Aj) 
we have \pi(Ai) — p, i (A - 2 )| < e(l + |/Aj(A 1 )|), where e = e 1 ^ — 1. 

Proof. Let A\ satisfy the hypotheses of the lemma, and let A 2 be a self-adjoint operator with 
||.Ai — A 2 H^i < S, where without loss of generality we assume that S < 1; define A(t) = A\ + 
t(A 2 — A\), for t £ [0,1]. Let A,(t) be the analytic branch of the eigenvalues of A(t) emanating 
from Ai(0) = pi(Ai), and denote by 4>i(t) a corresponding analytic eigenstate. By hypothesis 

||(^4 2 — Ai)<j>i(t)\\ < 5||A 1 0j(f)|| + 6 

< <5|Ai(t)| + <5t||(A 2 — A 1 ))^j(t)|| + 5 , 

which implies that || (A 2 — (|A*(t)| + 1). 

From 

|Ai(t)| = | (<f>i(t),(A 2 - Ai)<k(i))| < ||(A 2 - Ai)0i(t)|| 
and Gronwall Lemma we easily get 

|Ai( 0 )-Ai(l)| < (e^ -l)(|Ai( 0 )| + l), 

and then the thesis. □ 

When we consider parameterized families of self-adjoint operators, we can prove some prop¬ 
erties concerning the differentiability of spectral projections associated with separated portion of 
the spectrum. The statement here below deals with affine families, but can be generalized to more 
general settings (see e.g. nans] for similar arguments). 

Proposition A.3 Let K 0 be a self-adjoint operator, y be a Banach space with norm || • ||y and 
K(-) be a linear and continuous operator from y to the space of Kq- bounded self-adjoint operators, 
endowed with the norm || j|if 0 • Let moreover qo £ y and I C R. be an interval whose boundary points 
belong to the resolvent set of Kq + A'(go)- Then the spectral projection Pi(q) on I associated with 
the self-adjoint operator Kq + K(q) is well defined and (Frechet) differentiable on a neighborhood 
of qo- 

Proof. By assumption we have that ||A'(g) — K (go)|| k 0 < Cj|g — golly for some C > 0. Therefore 
for g in a sufficiently small neighborhood of go, C belonging to the resolvent set of K 0 + K(q 0 ) and 
setting R((q) = R(K 0 + K(q), £), we can write, thanks to ( |35| , 

R((q) = Rc(qo)(id+ (K(q) - K(q 0 ))R c (q 0 fj = A c (g 0 )^] ((AT(g 0 ) - K(q))R c (q 0 )) k . 

k—0 

Thus, from 

OO 

R d<l) - R d<lo) - A c (go)(AT(go) - A'(g))A c (g 0 ) = A c (g 0 )^ ((A'(g 0 ) - K(q))R c (q 0 )) k 

k=2 


we conclude that there exists a constant C > 0, continuously depending on such that 

\\Rdq) - Rd<lo) - R((qo)(K(q 0 ) - A'(g))A c (go)|| < Cj|g - g 0 ||y, 
which guarantees that 

lim tt- 1 — | — 

ll«-9o||y->o ||g — g 0 ||y j r 


j> (Rdo) ~ R ddo) - R d<lo)(K(q 0 ) - K(q))Rdqo)jdC 


= 0 , 
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where T is closed curve in C enclosing I (and not containing any other element of <j(Kq + K(qo))). 

□ 


In the last part of this section, we focus on control-dependent Hamiltonians satisfying assump¬ 
tion (HO) and, when explicitly said, assumption (HI) too. 

First of all, it is easy to see that for any Ui, U 2 the norms II •HtfOn), II ' |k( U2 ) are equivalent, 
thanks to the Uo-smallness of the control Hamiltonians. 

Let us now focus on the eigenvalues of H( u). From Lemma [5j and the equivalence of the 
norms || • ||#( u ), the eigenvalues A^(-) of H(-) are locally Lipschitz, and the corresponding Lipschitz 
constants locally depend on the magnitude of A,;(-). 

Let u be a conical intersection between the eigenvalues A j and Aj+i, that satisfy a gap condition, 
according to (HI). By the definition of conical intersection and the Lipschitz continuity of the 
eigenvalues we can conclude that there exist a suitably small neighborhood U of u and two 
constants C\ > 0 and C 2 > 0 such that 


Aj+i(u) — Aj(u) > Ci|u — u| V u G U 


(37) 


and 


|Aj(u) — Aj(u 7 )| < C 2 |u — u'| V u,u' e U, i = j,j + 1. (38) 

Moreover, if we consider two eigenvalues Aj, Aj+i, possibly intersecting, and isolated from the 
rest of the spectrum, the projection P u associated with these two levels is smooth with respect to 
u. The result holds also for any portion of the spectrum of H( u), in presence of a gap (see [3D]). 

On the other hand, the projections Pj , Pj+i, associated respectively with A j and Aj+i, are 
smooth with respect to u outside the singularity, while the presence of the conical intersection 
determines a lack of continuity at u. Nevertheless, along regular curves passing through the 
singularity, it is possible to extend these projections, obtaining operators whose regularity depends 
on the one of the curve, as stated in the following result. 

Proposition A.4 Let 7 : / —>• K 3 , I = [— R, 0], be a C k (I) curve such that 7 ( 0 ) = u is a conical 
intersection between the eigenvalues A j and Xj+i and 7 (f) 7 ^ 0 for every t £ I, and consider its 

k-jet at the origin £k(t) = 7 ( 0 ) + ^ 7=1 jk ^r 7 (t)|t=o- Then Pji'yi-)) is C k on [—R,0), it is C k 
at the singularity, and 

lim Wi= Dm 4f p i(4(i)), l = 0,..., k - 1, 

t—¥0~ dt l t-¥ 0 - at 1 

where the limit above holds in the operator norm. The same result holds for Pj. pi(7(0) - 


Proof. We first consider the case k = 1. Without loss of generality we assume I'V (0) | = 1. Let 


p= C i/4, where C\ is as in (37), and for every t £ [—R, 0) consider the circle T t C C of radius pt 


centered at Xj(j(t)). For a set A C C, we denote by d(z , A) = inf^g^t \z — x\ the distance between 
the point z and the set A. There exists 0 < T < R such that for every t £ [—T, 0) 

lk+i(7k)) - Aj(7(t))| > = 3 pt 


so that |Aj+i( 7 (i)) — Cl > 2 pt for every f £ T t . Thus d(£, cr(7J(7(t)))) = pt and, by (38) and 
the definition of fi(-), d((,a(H(£i(t)))) > pt/2, up to reducing T. Therefore from the classical 
identity holding for self-adjoint operators ||(X — £id) -1 || = d((, cr(A)) -1 (see e.g. JH]), for ( £ Tj 
there hold 


\R( 7(t),C)ll = 


1 


4C,cr(iL(7 (t)))) 


1 


= -, \\R{h{t),0\\ = 


1 2 

d(C,<r(H(£i(t)))) ~ pt ( ' 39 ^ 
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In order to get the thesis, we prove that 


lim 

t—y 0 


jf (Rh(t),0~R(h(t),0) dC 


= o. 


Estimate (35) gives 


Oil <c 


4 (*) ~l{t) 


for some C > 0, which, together with (34)-(39) and the definition of yields the thesis. 

Let us now tackle the general case; the proof follows similar arguments. We define the circuit 
r T as above, and we notice that for every fixed t £ (— T, 0) there is a neighborhood I T of r such 
that (39) can be replaced by the similar estimate 


\Rh(t),C)\\<- 

pt 


P(4(*),C)II < -i, 

pt 


(40) 


holding for every ( £ T r and t £ I T . For every t £ I T we have that 
pi} r r 


^T-PjW*)) = -(2wi) dC = -(2ni) 1 £ ^R^t), C) d,C, l<k- 1 


and thus, by applying (34), 

^p. 

dt l 3 


= -m- 1 jf o(H(£ k (t))-jT(7(*)))ij(7(t),o) 


dC: l < k—1. 


The proof can then be easily completed by applying recursively the identity 


the estimates (40) and (351, and by exploiting the regularity of 7 (-) and the definition of □ 


B Genericity of conical intersections 

In this section we discuss the occurrence of conical intersection for certain classes of Hamiltonians. 
More precisely, we show that conical intersection are not a pathological phenomena: indeed, in 
the finite dimensional case and for some specific families of infinite dimensional control-affine 
Hamiltonians, the fact of being a conical intersection is a generic property of eigenvalue crossings, 
with respect to the controlled Hamiltonians 0 

To prove this genericity property, in this paper we use transversality theorems that rely on 
the second countability of the family of Hamiltonian operators under consideration. In the finite¬ 
dimensional case, all these hypotheses are fulfilled. Concerning the general infinite-dimensional 
case, where we take the controlled Hamiltonians Hi, H 2 , H 3 as self-adjoint operators on H, clas¬ 
sical transversality theorems do not apply, since the space of self-adjoint operators is not second 
countable (even if we restrict our attention to bounded self-adjoint operators). However physi¬ 
cally relevant Hamiltonians often belong to particular families of operators which happen to be 
second-countable: we will focus on one of these cases. 

In the following, we will consider the class J- of Hamiltonians of the form (q) = I\q + K{q), 
where K$ is a fixed self-adjoint operator, q belongs to a Banach space y and K(-) is an injective 
linear and continuous operator from y to the space of i^Q-small self-adjoint operators, endowed 

1 We recall that a property is said to hold generically in a Baire space X if it is satisfied for all elements belonging 
to a residual subset of X, that is a set containing an intersection of (at most) countably many open and dense 
subsets of X. 
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with the norm || • ||x 0 . We will also assume that the parameterized family T satisfies the following 
condition called Second Strong Arnold Hypothesis. 

Second Strong Arnold Hypothesis (SAH2) : Assume that A is an eigenvalue of Jh(qo) for 
some qo £ y of multiplicity greater or equal than two. Then there exist two orthonormal eigenstates 
if}i,ip 2 °f A(f{qo) pertaining to A such that the three linear functionals 


hi - hi ■■ q (4>i,K(q)4’i) - ( 02 , K(q)ip 2 ) 
f 12 ■ q 1 t (0i, K(q)ip 2 ) 
hi ■ q < 02 , K{q)ipi) 


are linearly independent. Equivalently, the linear map $(•) = (/n(-) — /22(0, 9^ e (/i2( - )), 3'n(/i2( i ))) 
is surjective from y to R 3 . 

We call V the subset of y such that the Hamiltonians in Jh(V) have double eigenvalues. For 
every interval I and every open set IT in y , we denote by V 1 ^ the subset of elements in U such 
that the corresponding Hamiltonians have an eigenvalue in / of multiplicity two, isolated from the 
rest of the spectrum. 

In particular, under some additional regularity assumptions on the spectrum of the operators, 
(SAH2) guarantees that T> has codimension 3 in y. More precisely, for a sufficiently small interval 
I and a sufficiently small open set U the set V I,u is a smooth manifold of codimension 3 (see pUl). 

The conicity of eigenvalue intersections correspond to a geometric property in the space of 
parameters, as the following result shows. 

Lemma 6 Let H( u) = H 0 + U\Hi + U 2 H 2 + U 3 H 3 belong to T for every u £ R , and assume 
that it has an isolated double eigenvalue A at u = u; then there exists a unique q £ y such that 
H (u) = Jh(q) and unique qi,q 2 ,q 3 £ y such that Hi = K(qi), i = 1,2,3. Assume moreover that 
there is a neighborhood U of q iny and an interval I containing A such thatADis a submanifold 
of codimension three in y. Then u is a conical intersection for H(-) if and only if for every 
direction v £ R 3 the vector q v = V\q\ + v 2 q 2 + ^393 is not tangent to V I,U at q, that is, the affine 
space {q + q v : v £ R 3 } is transversal to V I U at q. 


Proof. The existence of q and qt as in the thesis comes directly from linearity and injectivity of 
K. 

If there exists some v £ R 3 such that q v is tangent to V I - u at q , then u cannot be a conical 
intersection, since in that case the distance between the eigenvalues intersecting at u is of order 
oft) along the line 1 1 —> u + tv. 

Let us now prove the converse statement. Denote by Ai(u) and A 2 (u) the two eigenvalues of 
H(u) crossing at u, with Ai(u) = A 2 (u) = A. 

Under the assumptions of the lemma, we can deduce the following facts. 


• Possibly reducing / (still containing A in its interior) and the neighborhood U of q, Jh(q) 
contains exactly two eigenvalues in /, counted with their multiplicity, for every q £ U. 


• Denoting with M(q) the sum of the eigenspaces of Ah(q) associated with the eigenvalues in 
I and with Pi{q) the orthogonal projection on M(q), we have that, possibly reducing U, 


\\P I (q)-P\\<l Vg £ IT , 


where P = Pi(q), and moreover Pi(q) is a differentiable function of q in IT, by Proposition A.l 
and Proposition |A.3| 




_ _ _ _I /n _ 

The map S(q) = Pi (g) (id + P(Pj(q) — P)P) P is an isometric transformation from M(q) 

onto M(q) (see e.g. [2H Section 105]), and it is differentiable with respect to its argument. 
Therefore the map 


f(q) = S{q) 1 Jh(q)S{q) 


is a differentiable mapping from U to the space of self-adjoint operators on M(q), and the 
eigenvalues of f(q) are the same as the eigenvalues of Ah ( q ) in I. 
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It is easy to see that V 1 - 11 C / _1 ({^icl : /x £ R}), where id denotes the identity on M(q). 

Let us now assume that the intersection between the eigenvalues Ai and A2 is not conical, that 
is there exists a unit vector vfl 3 such that 

A 2 (u + tv) — Ai (u + tv) = o(t), 

and we consider the curve N(t) = f(q + tq v ) in the space of self-adjoint operators on M(q), that 
we write as two dimensional Hermitian matrices in a basis made of eigenstates relative to Ai and 
A 2 ; it holds 



for some complex functions a(-), 6(-), c(-) satisfying a(0) = c(0) = A and 6(0) = 0. Since the eigen¬ 
values of N(t) coincide with those of if (u + fv) contained in J, it holds y^(a^ty^c(t))2~+4j6(t)p^ = 
o(t) and in particular, by the analiticity of Ai(u + tv) and A 2 (u + tv) with respect to t, it is easy 
to conclude that lim 4 _ >0 + (a(t) — c(t)) = 0 and lim t _ >0 + b(t) = 0, that is, N(t) is tangent to the 
space {/x.id : p, £ R} at the point Aid. 

Since by definition the family T satisfies the condition (SAH2), the map / is transversal to 
{/x id : /x £ M}. Then we can conclude that 

(D q f)-\{pid:pGR})=T q D I ’ U 

(see e.g. [TJ Corollary 17.2]) and, in particular, q v is tangent to D I ' U at q. □ 

B.l Finite-dimensional case 

The class of Hamiltonians under consideration is here J- = iu(n), that is the set of Hermitian 
n x n matrices. Trivially, in this case y coincides with T. It is well known (see e.g. [3j 1331 1 
that for sufficiently small / and U the set V I ' U is a smooth manifold of codimension 3 in ?'u(n), 
and, moreover, that the subset of Hermitian n x n matrices having at least one eigenvalue of 
multiplicity 3 has codimension 8 in iu(n). By second countability we can extract a countable 
family of pairs {{IkMk)}ke N that satisfy the property above in such a way that V = UkV Ik ’ Uk . 

Lemma 7 Fix H 0 £ iu(n). Borrowing notation from w, let us define the map p : iu(n) 3 —> 
C(R 3 \ {0}, iu(n)) as 

p(H u H 2 ,H 3 ) = H{-) (41) 

with H( u) = H 0 + U1H1 + u 2 H 2 + u 3 H 3 , and let ev p : (R 3 \ {0}) x iu(n ) 3 —?• iu(n) be defined as 
ev p (-, Hi, H 2 , H 3 ) = p(Hi, H 2l H 3 ). Then ev p is transversal to T> Ik,Uk for every k £ N. 

Proof. If ev p (u, H±, H 2l H 3 ) T> Ik ’ Uk , then the thesis trivially holds. Assume then that ff( u) = 
ev p (u, H ll H 2 , H 3 ) £ Uk has a double eigenvalue A £ Ik for u^0. 

The differential of ev p at (u, Hi, H 2 , H 3 ) along the direction (Su,SHi,SH 2 ,SH 3 ) is given by 

3 

Dev p\^,Hi,H 2 ,H 3 )^ u ’ SHl ^ H2 ' SH ^ = {uiSHi + SuiHi) . 

i=i 

Let us consider the three directions Vi = (0, 0, 0, uiai,u 2 ai,u 3 cri), where the operators 07, l = 1, 2, 3, 
are given by 

o-i = {1P2, -)vi + (V’l, -)V 2 cr 2 = i(tp 1, -)ip 2 - i(ip 2 , -)ipi 0-3 = (y>i, -)ipi - {ip 2 ,-)<P 2 , 

and ip 1 and ip 2 define an orthonormal basis of the eigenspace of H( u) relative to A. 

Let us consider the eigenvalues of 


3 

H{ u) + e Dev p \ (a Hi H2 Hs) [oiui + a 2 v 2 + a 3 v 3 ] = H{ u) + e|u| 2 ^ a^. 

i= 1 
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It is easy to check that the degenerate eigenvalues split and their difference is equal to 2e|u| 2 |aj. 
Therefore, span{Dev p | (Q [• vi\ : Z = 1,2,3} is a three dimensional space having trivial inter¬ 
section with T H ( a ' ) T> Ik ’ Uk . Since the codimension of V Ik ’ Uk is 3, this means that span{Dev p | ( _ H ^ [vt] : 

i = 1,2,3} is transversal to T m a) V Ik ’ Uk , and by definition of transversality of a map we get the 
proof. □ 

Proposition B.l Let H$ £ iu(n). Then generically with respect to (Hi, H 2 , H 3 ) £ iu(n) 3 , all 
double eigenvalues of H Q + u %Hi correspond to conical intersections. 

Proof. Define H( u) = Hq + X^=i u iHi- Thanks to Lemma [j] we can apply the Transversal 
Density Theorem !lj with W = D Ik,Uk for k £ N, and p(Hi, H 2 , H 3 ) = H(-), and obtain that the 
set of Hamiltonians (H\,H 2 , H 3 ) such that H(-), as a function defined in R 3 \ {0}, is transversal 
to V Ik,Uk is residual in i u(?r). In particular by applying Lemma [6] and taking the intersection of 
the previous residual sets over all k we get that, generically, all double eigenvalues of H( u) with 
u /0 correspond to conical intersections. 

Let us now consider the case where Hq has a double eigenvalue A, and let ifi and V ’2 define an 
orthonormal basis of the eigenspace relative to A. Consider the real-valued multi-linear map 

/^((Vq,#!^)) 3m«V>i,i?iV ; 2)) ((^ 2 ,# 1 ^ 2 } - (■0i,^iV’i))\ 
F(H 1 ,H 2 ,H 3 ) = det fHe((V’i,i? 2 'i/’ 2 )) 3m((ipi, H 2 if 2 )) (<^ 2 , # 2 ^ 2 ) - (ipi, # 2 ^ 1 )) . 

yffle ((4>i, H 3 'i/j 2 )) 3m((^i,ff 3 ^2)) ((^ 2 , # 3 ^ 2 ) - (i>i, # 3 ^ 1 ))/ 

Recall that u = 0 is a conical intersection for Hq + u iHi if and only if F(H\, H 2 l H^) 7 ^ 0. 
Moreover, being F continuous, we obtain that F~ 1 (M. \ {0}) is an open subset of iu(n) 3 . This 
subset is non-empty because the map 

H £ iu(n) >->• (fte ((ipu Hip 2 )), 3m ((%pi,Hif 2 )), ((^ 2 , Htp 2 ) ~ (ipi, Hipf)) ^ e R 3 

is surjective. The density comes directly from multi-linearity. □ 

B.2 Infinite dimension: the case of electromagnetic Hamiltonians 

Let us consider the class T of Hamiltonians of the form H = -A + L-i(VA + AV), where A 
denotes the Diriclilet Laplacian on a given Lipschitz bounded domain D C R 3 , V is a scalar con¬ 
tinuous real-valued function on its closure Cl, that should be thought as a multiplication operator, 
and A is a C 1 vector-valued real function from Cl to R 3 . We focus on this family of Hamiltonians, 
since they happen to be largely used to model quantum systems driven by electromagnetic fields. 

The self-adjoint operator Ha = -i(VA + AV) acts on the elements of its domain as follows: 

—i (VA + AX7)ip = — i A • Vip — i div(A')/)). 

Since, as it can be easily seen, the map (V, A) 1 —> V + Ha is linear and injective, T has a Banach 
manifold structure modeled on the space y = C(D,R) x C 1 (D,R 3 ); moreover, since Cl is bounded, 
it turns out that both y and F are separable, and thus second countable. 

It is not difficult to show, by a direct integration by parts and applying the inequality (^ 1 , ip 2 ) < 
^(£|| , 0i || 2 + jll^H 2 ), that self-adjoint operators of the form Ha are A-small. Therefore they can 
play the role of control Hamiltonians in our setting. Similarly it can be shown that each Ha is 
form-bounded with respect to — A (as a quadratic form, see E3 Chapter X]) with a relative bound 
that can be chosen smaller than one. Thus [25J Theorem XIII. 68 ] ensures that the Hamiltonians 
of the form —A + V — *(VA + AV) have compact resolvent so that their spectrum is purely 
discrete with a finite number of eigenvalues in each compact subset of R. Notice moreover that 
the topology induced by the norm || • ||a on F is coarser than the topology inherited from y. 

For the family of Hamiltonians F defined above we essentially repeat the same argument as in 
the finite-dimensional case to show a genericity property of conical intersections. Before stating 
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the main results of this section, some important remarks are in order. First of all we observe 
that the operator Hj± plays a crucial role for the existence of conical intersections for controlled 
Hamiltonian operators of the form Hq + u\H\ + U2H2 + u 3 H 3 belonging to J- . Indeed if one 
considers controlled operators belonging to the class of Schrodinger operators of the form —A + H 
with V a real-valued function, then conical intersections are never present. This can be seen as a 
consequence of the fact that the terms with i = 1,2,3, computed with respect to an 

appropriately chosen orthonormal basis of eigenfunctions {VbljeN °f — A + V, are real and thus 
the first two columns of each conicity matrix are equal. On the other hand, examples of conical 
intersections for controlled Hamiltonian belonging to the family T are not difficult to find, as 
shown below. 

Example. Consider the Hamiltonian H( u) = — A + UiVi +U2V2 + u 3 Hj±, where 

Vi(x) =xl + £3, V 2 (x) = x 2 x 3 , A = ( 0 , -X3/2, x 2 /2) t , 
x = (xi,x 2l x 3 ) e H = (0,1) x (0, a/ 3) x (0, a/5) 


with Dirichlet boundary conditions. We claim that H{ 0) (representing the potential well in H) 
admits conical intersections of eigenvalues. Indeed eigenvalues and eigenfunctions of H( 0) take 
the form 


A 


‘311 32,33 


(,? + f + f), 




n , 32 ,33 


, x 2 a /2 (ji'KX 2 \ . fj 3 nx 3 \ 

(x) = ^ smOr^r) sm ) sin ) 


where ji,j2,j 3 are strictly positive integers. Then one has that for instance Aip^ = Ai j2 ,2 corre¬ 
sponds to a double eigenvalue. A direct computation shows that the associated conicity matrix is 
nonsingular. 


Here we focus on the three-dimensional case since, unlike the other cases, it has a clear physical 
interest: in that case the vector A, called vector potential, is related to the action of a magnetic 
field B on the physical system determined by the relation B = V x A. However, let us observe 
that the fact that the domain Q is assumed to be a subset of R 3 is not crucial for the mathematical 
formulation of the problem and the correctness of the following results (the previous example, for 
instance, can be directly recast in a two dimensional setting since the variable x\ does not play any 
role there), unless the dimension is one. Indeed, in the latter case, it is easy to see that for every 
continuously differentiable A, the energy levels of the Hamiltonian — d 2 + V(x) + H\ coincide 
with those of the Hamiltonian —d% + V{x) — A 2 . It can be easily shown that Hamiltonians of the 
latter form, on bounded intervals and with Dirichlet boundary condition, do not admit degenerate 
eigenvalues. Therefore the results below do not provide any information in the one dimensional 
case. Note that the fact that for one dimensional systems a magnetic field can always be reabsorbed 
by an electric field is well known in physics. 

Let us now proceed with the study of the genericity properties of conical intersections for the 
class T of controlled Hamiltonians under consideration. 

First of all, we notice that the class T fits the formulation given at the beginning of this 
appendix, with .X(V. A) = —A + V + Hj±. We claim that the set V C y which parametrizes 
the operators admitting double eigenvalues has codimension three in y. The claim is proved once 
shown that all elements in T satisfy (SAH2) (see pH]). 


Lemma 8 Every element of J- satisfies (SAH2) for any multiple eigenvalue. Moreover, the 
restriction ( I > |{(y ! o):yeC(n,R)} (defined in the statement of ( SAH2)J has rank at least two. 

Proof. Let us consider K £ J, and assume that A is a multiple eigenvalue of H. By contradiction, 
assume that there exist three complex scalars a, b , c and two eigenstates of H relative to the 
eigenvalue A such that the functional 


a(/ n (V, A) - / 2 2 (H, A)) + bf 12 (V, A) + c/ 21 (V, A) (42) 
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is identically zero. Notice that this fact does not depend on the particular choice of the orthonormal 
eigenstates ip\, ip 2 - 

Integrating by parts the terms of the kind (ipi, taking into account boundary conditions 

on the eigenfunctions, we can write the functional above as 


VB - iA ■ £, 


where 


B = a|'0i | 2 - a|' 02| 2 + W'i '02 + cip\ip 2 , 

£ = a(^W>i - ipiVipl - ip 2 Vip 2 + V’sVV’s) + 6 (^i*V ^ 2 - foVipl) + c(^VV>i - ipiVip* 2 ). 


By arbitrariness of V and A, the expression ( |42[ ) is identically zero only if B and £ are identically 
zero on f 2 . 

Let us start by assuming that there exist two orthonormal eigenfunctions ip \, ip 2 such that 
\ipi\ = \ip 2 \ on £1 and three scalars a, &, c such that the functional (42) is zero. Without loss of 
generality we can assume that a £ R. In particular, from B = 0 we obtain that either ip 1 and 
ip 2 differ only by a constant phase, which contradicts their linear independence, or b = c = 0. In 
the latter case, denoting ipj = ipe l6j , j = 1,2, it turns out that £ = 2ia\ip\ 2 S7 (61 — 0 2 ). Then 
V(#i — 0 2 ) = 0 wherever \if>\ 0, and in particular 0\ — 0 2 is constant on a open set. Up to a 

phase change of the eigenfunctions, ip± — ip 2 is an eigenfunction which is null on a open set, which 
implies, by the unique continuation property (see e.g |19j). that ipi = ip 2l which is a contradiction. 
We can then conclude that there are no orthonormal eigenfunctions with equal absolute value that 


make (42) identically zero. 


Let us now assume the general case in which there exist ip 1 , ip 2 and a, b , c with a £ R such that 


the functional (42) is zero. Again from B = 0, we obtain that 


3m (B) = 3m(bip\ip 2 + cV'iV’S) = (b — c*)ip\ip 2 + (c — b*)ip\ip 2 = 0 a.e. on Cl. 


If b = c*, this condition is automatically satisfied and we write B as 

a|'0i| 2 - a\ij) 2 \ 2 + biplip 2 +cipxip 2 = (^*>$ 2 ) (l ■ 

We can then diagonalize this quadratic form, ending up with two orthogonal eigenfunctions <pi , (p 2 
of H, associated with A, such that 

v / a 2 T|feF(M 2 -1^ 2 | 2 ) = 0, 


which is not possible thanks to the arguments above. 

Let now b ^ c*. By unique continuation property, we can assume that ip*ip 2 is not identically 
zero, that implies that 

b — c* ipxip 2 
b* - c ip*ip 2 ’ 

which proves that the difference between the phases of ip 1 and ip 2 is constant on f 1 and, in partic¬ 
ular, it can be set to zero. This in particular leads to 3m(6) = —3m(c). Let us then set b = /3 + ir, 
c = 7 — ir, for some /3, 7 , r £ R. 

Let us write ipi = <pie l< ^ and ip 2 = <t> 2 ^■, for some real-valued functions <f> 1 , <p 2 and £. Then by 
computations it follows from £ = 0 that 

03 - 7)(0iV& - fcV 0 i) + 2 i(r(^ 1 V& - foVfa) + a(VC)(</>? - <pl) + (7 + P)(VQ4>i<h) = 0. 

By direct computation one checks that <p iV0 2 — <p 2 ^<Pi is proportional to ipiVip 2 — V?/q, and 
applying again the unique continuation property it turns out that the latter cannot be identically 
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zero on an open set in Q. Therefore it must hold /3 = 7 , that is c = b*. By contradiction, we see 
that (SAH2) is verified. 

The proof of the second statement follows the same arguments and is thus omitted. □ 

We remark that, thanks to the fact that the spectrum of any operator in T is discrete with 
no finite accumulation points, and that the eigenvalues are continuous with respect to the pair 
(V, A) (see Lemma [HJ, then for every (V, A) £ y such that J^(V, A) has a double eigenvalue A 
there exist a neighborhood I of A and a neighborhood U of ( V , A) such that the subset V I U is a 
smooth submanifold of codimension three in y. In particular, as in the finite dimensional case we 
can find a countable family {(IkMk)}k such that T> = UkT >Ik,Uk 1 with V Ik ’ Uk smooth submanifold 
of y of codimension three. 

Let us now consider controlled Hamiltonians in T of the kind H 0 + U\X\ +U 2 V 2 + u 3 Ha, where 
H 0 £ F, 14,14 £ C(fi,R) and A £ 

Lemma 9 Let H = — A + fo + Ha, for some V £ C(f2,R) and A £ C 1 (f2,R 3 ). Let X = {u £ R 3 : 
m(+« 2^0 and u 3 7 ^ 0} and let us define the following m.ap p : (C(fi,R )) 2 x C 1 (H,R 3 ) —> C(X, y) 
as 

p(V 1 ,V 2 , A)(u) = (mVi + u 2 V 2 , u 3 A). (43) 

Let ev p : X x (C(fi, R)) 2 x C 1 (H, R 3 ) —7 y be defined as ev p (-, 14, V 2 , A) = p(V 1 , V 2 , A). Then ev p 
is transversal to T> Ik,Uk for every k £ N. 

Proof. Let us fix some notations: we set K(V, A) = V + Ha and Jtf (V, A) = H + K(V, A). 

If ev p (u, Vi, V 2 , A) </ x> Ik ’ Uk , then the thesis trivially holds. 

Assume then that H( u) = Jff (ev p (u, V t . V - 2 . A)) has a double eigenvalue A £ Ik for some u £ X 
and (14,14, A) with ev p (u, Vi, V 2 , A) £ Uk, and that ip 1 , 1/12 are two orthonormal eigenstates of 
H( u) pertaining to A. Without loss of generality, we assume that Hi 7 ^ 0. Lemma [ 8 ] ensures that 
the rank of the map 

(V, A) ha ((^ 7 , (V + H A )ih), (V’i, (V + H A )ip2)*, ( 1 P 2 , (V + Ha)iPi) - (ipi, (V + Ha)iPi)^ 

is three, and its restriction to the space {(V, 0) : V £ C(f2,R)} has rank at least two. Then we can 
find three functions SU,5W £ C(f2,R) and S A £ C 1 (f2,R 3 ) such that the conicity matrix 

/ (ip! ,SUlp 2 ) (1 p!,6Ulp2 )* (lp2,SUlp 2 ) - (lp!,SUlp!) \ 

M(ip!,ip 2 ) = (ip 1} 6Wip 2 ) (ipi, 6Wip 2 )* {iP2,SWiP2)~(ip!,6WiP!) \ 

\{lpl,H S Alp2} (ipliHsAlfa)* (i>2,H5Alp2) ~ (4>1, Hs A^l )) 

is non-singular. In particular, this means that v = 0 is a conical intersection for the Hamiltonian 

if (v) = if (u) + ViSU + v 2 SW + v 3 H$ a- 

Moreover, we can write the Hamiltonian if(v) as 

3 

H(v) =H(u)+J2 v o K (Dey pka<Vl ,v 2 , A ) H'O 

3 =1 

3 

= H + K(uiVi + u 2 V 2 ,u 3 A) + v j K ( Dey P U a , Vl ,v 2 , a) K'O 

i =1 

where -Dev p | ( _ A) [«;] denotes the differential of ev p at (u, 14,14, A) evaluated on the variation 
w £ R 3 x (C(H,R )) 2 x C 1 (H,R 3 ), and the variations Wj are 

ivi = [0,0,0,(517/111,0,0], 

w 2 = [o,o,o,$wy« 1 ,o,o], 

w 3 = [0,0,0, 0 ,0, 8 X/u 3 ] . 
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Therefore, since the codimension of V Ik,Uk is three, and applying LemmajGjwith qi = Dev p | ( _ v2 A) [wi] 
and q = (u\Vi + U 2 V 2 , W 3 A), we get that 

^ = {« v :v€K 3 } + TqV Ik ,Uk C Im(De V(9 | (Q ^ V2 A) ) + TqV Ik ' Uk , 

and this concludes the proof. □ 


Proposition B.2 Let Hq G T and X as in the lemma above. Then generically with respect to 
(Vi, V 2 , A) G (C(f2,K)) 2 xC 1 (f2,K 3 ) ; all the double eigenvalues ofH(\T) = HQ+U 1 V 1 +U 2 V 2 +U 3 HA 
with u G X correspond to conical intersections. 


Proof. Thanks to Lemma [9] we can apply the Transversal Density Theorem and obtain that the 
set of triples (Vi, V 2 , A) such that p{V\, V 2 , A), as a function defined on X, is transversal to V Ik,Uk 
(at p(Vi, V 2 , A)) is residual in (C(Cl, K )) 2 x C 1 (fi, R 3 ). We can conclude as in Proposition B.l that, 


generically, all double eigenvalues of H( u) with u G X correspond to conical intersections. 
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